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from  NOAA  1TOS  satellites  by  means  of  non-linear  equations.  Equations  were 
also  developed  to  estimate  total  cloud  mass  and  thickness  from  the  same 
satellite  data. > The  visible  data  were  normalized  to  approximate  the  reflection 
geometry  of  both  sun  and  satellite  directly  over  the  cloud.  The  areal  mean 
temperature  from  the  infrared  imagery  data  was  found  to  be  the  most  signifi- 
cant single  predictor  for  erosion  parameters;  however,  the  most  accurate 
estimates  of  erosion  parameters  came  from  a combination  of  infrared  and 
visible  data.  Another  set  of  equations  was  developed  using  only  infrared  data 
and  climatological  temperature  altitude  profiles,  and  these  equations  were 
used  for  areas  without  sunlight  or  with  highly  reflective  snow  cover. 


^Extensive  applications  and  tests  were  made  of  the  equations  to  estimate 
erosion  parameters  from  satellite  data.^  One  year  (April  1973  to  March  1974) 
of  data  from  NOAA  polar -orbiting  satellites  was  processed  and  climatologies 
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radar  estimates  of  weather  erosion  and  also  a comparison  to  subjective  esti- 
mates of  erosion  based  on  standard  weather  reports. 


Unclassified 


>ICU"ITY  CLASSIFICATION  OF  THIS  FAOCnni.A  Dm „ tnimrmd) 


Preface 

The  following  report  relies  on  measurements  from  complicated  satellite  and 
aircraft  systems.  Technique  development  and  the  resulting  satellite  estimates  of 
reentry  vehicle  erosion  would  not  have  been  possible  without  the  assistance  of  many 
individuals  and  organizations.  Aircraft  sampling  of  clouds  underneath  satellite 
passes  was  directed  by  the  Convective  Cloud  Physics  Branch  (LYC),  AFGL,  under 
Dr.  Cunningham  and  later  Dr.  Barnes.  Captains  Brooks,  Cameron,  Delgago, 

Kelley  and  Lt.  Col.  Church  served  as  flight  directors  for  most  of  the  missions. 
Additional  flights  were  made  under  contract  by  Meteorology  Research,  Inc. , and 
Aeromet,  Inc.  Aircraft  data  reduction  was  provided  by  Mr.  Glass  of  LYC  and  by 
Digital  Programming  Services,  Inc.  Satellite  data  and  much  helpful  discussion  were 
provided  by  the  National  Environmental  Satellite  Service  and  National  Climatic 
Center  of  NOAA.  Helpful  studies  were  provided  by  Mr.  Cozzens  and  others  of 
Regis  College  to  display  estimates  of  cloud  properties  derived  from  satellite  data 
and  apply  the  estimation  techniques  to  geosynchronous  data  in  the  AFGL  MclDAS 
system.  Mr.  Myers  of  LY  and  the  MclDAS  team  directed  the  real  time  applications 
in  this  report.  Ms.  Bench,  formerly  of  the  Satellite  Meteorology  Branch,  wrote 
most  of  the  computer  programs  for  the  analysis  of  archived  satellite  data  with  assist- 
ance from  Envir  omental  Research  and  Technology,  Inc.  (ERT).  Radar  data  were 
taken  and  analyzed  by  the  Radar  Meteorology  Branch  under  Mr.  Glover  and  by  ERT. 
Important  consultation  was  provided  by  Mr.  Lund,  Design  Climatology,  and  by 
Dr.  Cunningham,  Senior  Scientist  in  Clo  Physics.  Finally,  every  step  of  the 
work  by  the  authors  was  encouraged  and  assisted  by  Mr.  Touart,  Manager  of  the 
AFGL  Weather  Erosion  Program. 


I 


3 


Contents 


1.  INTRODUCTION  9 

9 
10 
11 
11 


2.  OBSERVATIONS  14 

2.  1 Satellite  Data  14 

2.  1.  1 ITOS  Satellite  Data  15 

2.  1.  2 GOES  1 and  SMS- 2 Geosynchronous  Satellite  Data  23 

2.  2 Aircraft  Data  25 

2.  2.  1 C-130  Aircraft  and  Data  Reduction  25 

2.  2.  2 Lear  36  Jet  Aircraft  28 

2. 3 Radar  Data  28 

2.  3.  1 SPANDAR  Radar  28 

2.  3.  2 TRADEX  Radar  29 

2.4  Sample  Case— 8 October  1976  29 

2. 4.  1 Pre-flight  29 

2.  4.  2 General  Synoptic  Situation  30 

2.4.3  Hydrometeor  Mass  Profiles  32 

2.4.4  Satellite  Data  for  8 October  1976  37 

2.  5 Summary  of  Satellite  and  Erosion  Parameters  for  all  Cases  39 

3.  RELATIONSHIP  BETWEEN  SATELLITE  PARAMETERS  AND 

EROSION  PARAMETERS.  LIQUID  WATER  CONTENT  AND 

CLOUD  THICKNESS  42 

3.  1 Transformation  of  Satellite  Data  to  a Linear  Form  42 

3.  2 Dependent  Sample  Size  43 

3.  3 Erosion,  Liquid  Water  Content,  and  Cloud  Thickness 

Predictions  from  Brightness  and  IR  Data  44 

3.  4 Erosion,  Liquid  Water  Content,  and  Cloud  Thickness 

Predictions  from  IR  Data  Alone  46 

3. 5 Probability  Predictions  50 


5 


1.  1 The  Significance  of  Weather  Erosion 
1.  2 Previous  Erosion  Climatologies 
1.  3 The  Utility  of  Satellite  Data 

1.4  Related  Satellite  Studies 

1. 5 Approach 


Contents 


4.  APPLICATION-CLIMATOLOGICAL  AND  REAL  TIME 

EROSION  PARAMETERS  51 

4.  1 Program  Schedule  51 

4.  2 Grid  Size  5 1 

4.  3 Environmental  Severity  Index- Precipitation 

(ESI  ) Maps  52 

4.  4 ES1b  f<Pr  the  Kwajalein  and  Wallops  Test  Ranges  52 

4.  5 Realtime  Application  of  ESlp  Erosion  Parameter  52 

5.  VERIFICATION  AND  INTERCOMPARISON  56 

5. 1 General  56 

5.  2 Statistical  Approach  57 

5.  3 Radar  vs  Satellite  ESI„  Values  59 

5.  3.  1 SPANDAR  RADAR  vs  GOES-East 

Satellite— 24  February  1977  59 

5.  3.  2 TRADEX  RADAR  vs  GOES-West 

Satellite  — 8 July  1977  61 

5.4  Aircraft  vs  Satellite  ESln  Values  Near  Kwajalein  61 

5.4.  1 The  8 August  1976  Case  61 

5.4.2  The  12  August  1976  Case  65 

5.4.3  The  14  August  1976  Case  67 

5.  5 Effect  of  Slant  Viewing  from  a Satellite  69 

5.  6 Satellite  vs  AFGL-11  Erosion  Estimates  72 

5.  7 Conclusions  on  Verification  75 

REFERENCES  79 

APPENDIX  A:  Computer  Processing  of  Archived  Satellite 

Data  for  Climatological  Applications  83 


Illustrations 


1.  Satellite  Archive  Data  Array  Over  the  Northern  Hemisphere  17 

2.  Satellite  Archive  Data  Mesh  Size  vs  Latitude  17 

3.  Visible  Sensor  Imagery  Generated  from  the  Digital  Archive 

for  14  and  15  February  1975  18 

4.  Geometry  of  Reflectance  19 

5.  Pattern  of  Satellite  Data  by  Orbits  for  the  Northern  Hemisphere 

14  and  15  February  1975  21 

6.  Atmospheric  Probes  Aboard  the  C-130E  Aircraft  26 

7.  Precipitation  Hate— Mass  Density  Relationships  for  Rain, 

Small  Snow,  and  Large  Snow  27 

8.  Display  of  NOAA-5  Digital  Visible  Data  1512Z  8 October  1976, 

Surface  Synoptic  Pattern  at  1500Z  and  Location  of 
Reference  Points 


31 


Illustrations 


9. 

10. 

11. 

12. 

13. 


14. 


15. 

16. 
17. 


18. 


19. 


20. 


21. 


22. 

23. 

24. 


25. 


26. 

27. 

28. 

29. 

30. 


PPI  Display  of  Weather  Echoes  at  Satellite  Pass  Time 

Particles  Sampled  by  the  2D-PMS  Probe  During  Ascent 
8 October  197  6 

Vertical  Profiles  of  Mass  for  Particles  Greater  than  50  pm 
in  Diameter  During  Ascent  and  Descent  8 October  1976 

Vertical  Profiles  of  Mass  for  Particles  Less  than  50  pm 
in  Diameter  Measured  by  the  Johnson-Williams  Sensor 
During  Ascent  and  Descent  8 October  1976 

Character  of  the  Cloud  Tops  Observed  from  Top  Level  of 
10.  0 km  During  Ascent 

Thinning  Ice  Clouds  and  a Patch  of  Liquid  Cloud  Photographed 
from  12.  96  km  at  Start  of  Descent 

Aircraft  Sounding  Track  and  Area  of  Digital  Visible  Satellite 
Data  Assumed  to  Have  Passed  the  Track  Dui  ng  Ascent 
and  Descent 


Same  as  Figure  15  Except  for  IR  Data  in  °K 

Data  Points  and  Hyperbolic  Equation  Used  to  Relate  Visible 
Data  to  ESI 

P 

Plot  of  Equation  Used  to  Predict  ESlp  from  B^.  andlR,  Based 
on  4 1 Cases 

Plot  of  Final  Equation  Used  to  Predict  ESIp  from  Hj^ 

Probabilities  of  Exceeding  an  ESIp  of  5 Related  to  Satellite 
Measurements 

ESIp  in  the  Vicinity  of  Wallops  Island  1 April  1973,  from  the 
NOAA- 2 Satellite 

ESIp  in  the  Vicinity  of  Kwajalein  Atoll  9 April  1975  from  the 
NOAA -2  Satellite 


ESIp  (in  tenths)  in  the  Vicinity  of  Kwajalein  Atoll  21  July  1977 
from  the  GOES-West  Satellite 


ESIp  (in  tenths)  Derived  from  GOES-W 

Cumulonimbus  Anvil  Showing  Two  Levels  of  Spreading  Below 
the  Top 

Left:  Mass  Profiles,  at  the  Storm  Core,  3 and  6 km  NE  and  SW 
from  the  Core  Along  the  Aircraft  Track,  right:  Cross  Section 
of  the  Storm  Along  the  Aircraft  Track 

WSR-57  Radar  Echo  (stippled)  of  the  Storm  at  Satellite  Pass  Time, 
ESI  Values  Along  the  Aircraft  Track  and  Estimated  ESIp 
Comours  F 


Mass  Profile  and  Observer  Notes  Through  Dying  Cumulonimbus 
(Cb)  Near  Kwajalein  on  12  August  1976 

ESL  (in  tenths),  WSR-57  Radar  Echoes  (stippled)  and  Aircraft 
Sounding  Area  (heavy  box)  on  2143Z  12  August  1976 

Mass  Profile  and  Observer  Notes  Through  Active  Cb  Near 
Kwajalein,  2103-2142Z,  14  August  1976 


31  < 

33 

34 

35 

36 

37 

38 

39 

43 

46 

50 

50 

53 

54 

56 

62 

63 

63 

64 
66 

67 

68 


7 


Illustrations 


31. 

ESIp  (in  tenths),  WSR-57  Radar  Echoes  at  2131Z  (stippled)  and 
Aircraft  Sounding  Area  (heavy  box) 

69 

32. 

ESIp  (in  tenths).  Pattern  Calculated  from  GOES  W.  17° 

Zenith  Angle  to  Satellite 

71 

33. 

ESIp  (in  tenths)  Pattern  Calculated  from  GOES  E,  56° 

Zenith  Angle  to  Satellite 

7 2 

34. 

Method  1 Exceedance  Statistics  for  AFGL  11  and  Satellite  ESI 
V'alues  " 

74 

35. 

Method  2 Exceedance  Statistics  for  AFGL  11  and  Satellite  ESI 
Values 

75 

AI. 

Sequence  of  Computer  Programs  Used  to  Produce  Erosion 
Climatologies  and  Maps  from  Archived  Satellite  Data 

85 

Tables 

1. 

Conversion  of  1R  Coded  Values  to  °K  for  the  NOAA  Scanning 
Radiometer  Archive 

16 

2 

Correction  Factors  (\/ r)  to  Adjust  Observed  Brightnesses  to 
Reflectance  from  an  Isotropic  Surface  With  the  Sun  Overhead 

19 

3. 

Particle  Classification,  8 October  1976 

34 

4. 

Dependent  Data  Sample 

41 

5. 

Comparison  of  the  Correlation  Coefficients  (r)  and  the  Standard 
Errors  of  Estimates  for  ESI  Equations  Based  on  41  and 

42  Cases  ^ 

44 

6. 

Equations,  Correlation  Coefficients,  and  Standard  Errors  of 
Estimate  for  Prediction  of  Erosion  Parameters,  LWC  and 

Cloud  Thickness,  Based  on  41  Cases 

45 

7. 

Erosion  Parameter  Uncertainties  Corresponding  to  Maximum 
Monthly  Standard  Deviations  of  Temperature  at  Two  Altitudes 

47 

8.  Conversions  from  Equations  Using  Ih-pOnly  to  Equivalent  Equations 

of  TR  and  15^  49 

9.  Final  Equations  to  Predict  Erosion  Parameters  from  IL-p  Alone, 

Based  on  41  Cases  49 

10.  Program  Schedule  for  Satellite  Estimates  of  Erosion  Parameters  51 

11.  Frequency  of  Satellite  Coverage  at  Test  Kanges  During  One  Year  55 

12.  Root- Mean-Square  (rms)  Residuals  of  ESIp  in  Dependent  (40  case 

cample)  and  Independent  (2  X 20  case  samples)  58 

13.  Comparison  of  ESIp  Estimates  from  GOES- East  Satellite  and 

from  SPANDAR  Radar  60 

14.  ESIP  Determined  by  Aircraft  and  by  Satellite  in  Relation  to  the 

Oare  of  a Small  Isolated  Cumulonimbus  Cloud  65 


8 


L. 


J 


Estimates  From  Satellites  of  Weather  Erosion 
Parameters  for  Reentry  Systems 


1.  INTRODUCTION 

1 . 1 Hir  Si^nificancr  of  W ratlin  Eroaion 

The  nosetips  of  reentry  vehicles  ;ire  sub  jected  to  extremes  of  aerodynamic 
heating  and  pressure  while  decelerating  in  the  atmosphere.  The  air  temperature 
in  the  boundary  layer  often  reaches  values  above  10,000°F,  which,  combined  with 
the  surface  shear  forces,  promotes  intense  convective  heat  transfer  to  the  surface. 
Nosetips  made  of  material  with  high  heats  of  ablation  have  been  designed  to  provide 
protection  to  the  essential  parts  of  the  vehicle.  The  amount  of  nosetip  material 
required  for  heat  protection  during  clear  air  reentry  has  been  minimized  due  to 
cost  and  weight  considerations.  Additional  mass  toss  may  occur  due  to  hypervelo- 
city impact  with  atmospheric  hydrometeors.  Additional  nosetip  material  is  re- 
quired to  assure  vehicle  protection  margins  are  maintained.  Unpredictable  nose- 
tip mass  loss  may  also  affect  system  targeting  accuracy.  The  problem  of  hydro- 
meteor erosion  increases  with  velocity  and  must,  therefore.  Vie  considered  in  the 
design  of  high  performance  reentry  systems. 

Although  hydrometeor  erosion  is  only  one  constraint  among  many  in  the  design 
of  reentry  systems,  it  has  been  investigated  extensively  by  the  Advanced  Ballistic 
Reentry  Systems  (ABRF.SV  program  of  the  -\ir  Force  Space  and  Missiles  Systems 
Organization (SAMSOl.  Since  1973,  ABRES,  which  is  the  DoD  lead  agency  for  the 

(Received  for  publication  28  November  19771 
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design  of  reentry  systems,  has  sponsored  projects  to  develop  climatologies  of 
erosive  weather  over  target  areas.  This  report  summarizes  the  results  of  the 
Satellite  Correlation  Study,  which  uses  data  from  weather  satellites  to  estimate  the 
erosion  potential  of  clouds  over  target  areas.  Useful  by-products  of  this  study 
are  climatologies  and  real-time  estimates  of  erosion  potential  at  test  ranges. 


1.2  hrtiuui  Ihwiuii  (limaliilopn 

Climatologies  have  been  leve'oped  in  three  earlier  studies.  Theoretical  and 
empirical  methods  to  measure  cloud  and  precipitation  moisture  were  thoroughly 

reviewed  bv  Smith,  ' who  developed  relationships  to  specify  particle  water  or  ice 

2 

content  at  points  in  time  and  space.  Feddes  applied  these  relationships  to  the 
USAF  global  3DNEPH  analysis  output  to  produce  estimates  of  condensed  moisture 
content.  AFGL  experience  with  the  results  indicated  occasional  large  inconsisten- 
cies with  surface  observations  of  precipitation,  especially  in  the  cold  season.  In 
addition,  the  vertical  resolution  of  the  upper  layers  was  often  too  coarse  for 
hydrometeor  profiles  suitable  for  erosion  studies.  In  another  climatology,  known 
as  AFGL  1.  moisture  profiles  were  modeled  from  upper  air  (RAOB)  and  surface 
weather  reports.  The  model  was  applied  to  selected  RAOB  stations  in  Europe  and 
Asia.  Like  the  Smith-Feddes  model,  it  was  entirely  objective  and  computer  com- 
patible. A more  detailed  climatology,  known  as  AFGL  II,  was  developed  by  Cun- 

3 

ningham  and  Peirce  and  allowed  for  many  subjective  decisions.  Time-height  cross 
sections  were  prepared  for  11  stations  over  a one-year  period,  February  1973  to 
January  1974.  Preparation  of  the  cross  sections  is  described  by  Feteris,  Lisa, 

4 

and  Bussey.  Data  used  were  3-hourly  surface  reports  and  3DNEPH,  6-hourly 
RAOBS,  12-hourly  surface  maps  and  daily  DMSP  satellite  high  resolution  imagery. 
The  cross-sections  represented  a best  possible  cloud  depiction  prepared  by  experi- 
enced analysts.  Cloud  mass  densities  were  then  assigned  using  average  values  for 
cloud  types  and  temperature  modified  according  to  system  strengths  and  state  of 
development  or  decay.  Vertical  profiles  of  mass  density  were  recorded  every 
3 hours.  The  high  frequency  of  profiles  from  this  model  was  a distinct  asset  but 
its  practical  limitation  to  small  "areas  of  representativeness"  was  a weakness.  The 

1.  Smith,  R.C.,  Capt  (1974)  Atmospheric  Moisture  Parameterization,  ITSAFETAC 

TN74-1.  Washington,  D.C.  19  pp. 

2.  F’eddes,  R.  G.  , Capt  (1974)  A Synoptic -scale  Model  for  Simulating  Condensed 

Atmospheric  Moisture.  ItsAPfeTAC  TN  74-1,  Washington,  D.  C . , 21  pp 

3.  Cunningham,  R.  M.  , and  Peirce,  R.  M.  (1974)  Environmental  Definition  Plan 

Cross  Section  Analysis.  AFGL  (unpub.  ms.  ). 

4.  Feteris,  P.  J.  , Lisa,  A.S.,  and  Bussey,  A.  J.  ( 1975)  Environmental  Definition 

Program  Cross  Sectional  Analysis;  Summary  of  PatiTand  Analysis  Tech- 
niques,  \ Ec Hi- -TR -7 6-0002.  Final  llept.  F 19628 -74 -C -007  3,  Envlronme nt a 1 
lies,  and  Tech.  , Inc.,  Concord.  Massachusetts  01742,  34  pp. 
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limitations  were  addressed  by  Gringorten  who  developed  stochastic  models  ' 
approximating  the  occurrence  of  weather  events  various  distances  from  thr  t» 
of  the  event. 

1.3  The  Utility  of  Satellite  Data 

In  the  Satellite  Correlation  Study,  infrared  and  visible  satellite  data  hav.  • 
used  to  estimate  erosion  parameters  directly,  rather  than  the  vertical  profile 
hvdrometeor  mass  density.  The  erosion  parameters  are  weighted  integral- 
density  profiles  which  estimate  total  erosion  or  range  foreshortening.  1 
have  been  calculated  once  per  day  for  one  year  (April  1973  to  March  1974  1 .< 

of  Europe  and  Asia  and  for  test  ranges  at  Kwa.ialein  and  Wallops  Island. 

The  satellite  estimates  have  some  advantages  when  compared  to  the  ear  • 
studies.  First  of  all,  the  satellite  estimates  provide  areal  coverage  inforr  «'i  ■ 
so  that  it  is  known  when  widely  separated  targets  are  simultaneously  cover*  • • 
erosive  clouds.  Second,  the  estimates  do  not  rely  on  data  from  foreign  wear  • 
services  but  use  data  fr  m USA  satellites.  Third,  the  satellite  estimates  tr> 
sitive  to  the  upper  regions  of  clouds,  which  account  for  much  of  the  nosetu  • 
and  are  not  well  described  by  conventional  weather  data.  Finally,  the  sate! 
and  computational  techniques  lend  themselves  to  computer  processing  and  • • . 
time  applications.  Computer  processing  is  a distinct  advantage  in  dealing  w ' 
areas  of  Furope  and  Asia  once  per  day  for  one  year.  Real-time  application- 
useful  in  short  term  forecasts  for  the  Wallops  Island  Range  and  the  Kwaiale  '» 

sile  Range,  since  there  are  few  reports  of  cloud  systems  over  water.  I 

1.4  Related  Satellite  Sludie* 

Recent  publications  in  satellite  meteorology  and  in  the  radiative  properties 
clouds  have  been  reviewed  for  sources  of  input  to  weather  erosion  requirer  < •- 
The  requirements  are  severe.  The  most  useful  information  for  nosetip  from 
models  is  the  vertical  profile  of  hydrometeor  mass  density.  It  is  extren  civ  u - • 

Iv,  however,  that  detailed  vertical  profiles  could  be  inferred  from  present  -.it* 
measurements  since  present  sensors  are  passive  instruments  recording  r«  fie.  (• 
sunlight  or  thermal  emission  from  clouds.  The  sensors  are  sensitive  to  the 
cal  depths  of  clouds,  which  do  not  bear  a one-to-one  correspondence  with  vert 
profiles  of  cloud  mass  since  many  different  profiles  could  have  the  same  opt 
depth.  Consequently,  it  is  not  surprising  that  there  have  been  no  studies  to  r> 
trieve  detailed  vertical  profiles  of  cloud  mass. 


5.  Gringorten,  1.  1.  (197  6)  Areal  Coverage  Estimates  by  Stochastic  Modeling. 
AFGL -TR -7 6-0 148,  ERE  rto.  67  3.  bb  pp. 
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Fortunately,  estimates  of  integrated  profiles  of  hydrometeor  mass  density  are 
also  useful  for  nosetip  erosion  studies  and  satellite  measurements  are  useful  for 

g 

estimating  integrated  profiles.  In  research  sponsored  by  AFG1,,  I.iou  and  Stoffel 
7 

and  Feddes  ami  I.iou  have  reported  on  techniques  to  estiniate  the  total  mass  of 
cirrus  clouds  and  middle  level  clouds  from  infrared  data.  Estimates  of  integrated 

g 

liquid  water  from  microwave  data  have  been  reported  by  Staelin  et  al  and  r ainfall 
rates  by  Wilheit  et  al.  ' The  estimates  were  made  only  over  oceans  since  micro- 
wave emissivities  are  highly  variable  over  land. 

Although  these  techniques  provide  useful  information  for  some  cloud  conditions, 
they  are  not  presently  suitable  for  a climatology  of  nose  tip  erosion.  The  micro- 
wave  technique  is  a poor  predictor  of  ice  above  the  melting  level  of  a rain  cloud, 

but.  erosion  is  quite  sensitive  to  ice  particles  at  high  altitudes.  The  infrared  tech- 

6 7 

niques  developed  by  I.iou  ’ and  others  show  promise  for  thin  ice  clouds  and  some 
layered  situations  but  are  of  limited  application  for  thick  cloud  systems. 

These  techniques  are  based  on  models  of  cloud  radiation.  However,  aircraft 
observations  of  cloud  particles  in  heavy  weather  situations  hazardous  to  nosetips 
have  found  distributions  of  cloud  particles  which  are  extremely  unlike  distributions 
used  in  any  cloud  radiation  models.  The  aircraft  sensors  have  detected  ice  frag- 
ments, aggregates,  graupel,  rimed  particles,  dendrites,  and  other  irregularly 
shaped  particles.  These  uncertainties  as  to  shape  make  it  difficult  to  calculate 
phase  functions  for  single  scattering,  leaving  aside  the  enormous  additional  com- 
plications of  multiple  scattering.  The  aircraft  have  also  observed  great  variability 
of  cloud  particle  distributions  in  different  volumes  of  the  same  cloud,  which  implies 
that  a theoretical  approach  would  have  to  guess  at  the  particle  size  distribution  of 
the  cloud  and  consider  major  variations  with  altitude. 

The  difficulties  in  generalizing  theoretical  approaches  to  the  full  range  of  cloud 
conditions  erosive  to  nose  tips  led  us  to  review  the  literature  on  empirical  ap- 
proaches. We  found  numerous  publications  describing  methods  to  infer  cloud  pro- 
perties from  satellites.  Cloud  properties  which  have  been  estimated  from  satellite 
data  and  which  are  related  to  the  density  profile  of  hydrometeors  include  cloud  top 

6.  I.iou,  K-N,  and  Stoffel,  T.  (1976)  Remote  Sensing  of  Cirrus  Cloud  Composi- 

tions  From  Satellites,  A FGL-TR  -76-0027 , 81  pp.  University  or  Iltah, 

Sale  I.ake  City,  Utah  64112. 

7.  Feddes,  R.  G.  , and  I.iou,  K-N  ( 1977 ) Cloud  Composition  Determination  by 

Satellite  Sensing  Using  the  Nimbus  Vf  lligh  Resolution  Infrared  Sounder, 
AFC,L-Tn^7-im!r;  lAOpp,  TrnlvorsTfy  (TrntalT.  Salt  TUKc  niy,  Utah  844 12. 

8.  Staelin,  1).  H.  , Kunzi,  K.  F.  , Pettyjohn,  R.  I,.  , Poon,  R.K.L.  , and  Wilcox,  R.W. 

(197  6)  Remote  sensing  of  atmospheric  water  vapor  and  liquid  water  with  the 
Nimbus  5 microwave  spectrometer,  .).  Appl,  Meteor.  15: 1204-1214. 

9.  Wilheit,  T.T.  , Chang,  A.T.C.,  Rao,  M.S.  V.,  Rodgers.  E.  B. , and  Theon.  J.S. 

(1977)  A satellite  technique  for  quantitatively  mapping  rainfall  rates  over 
the  oceans,  .1.  Appl.  Meteor.  16:551-560 
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altitude,  cloud  type,  cloud  thickness  and  rainfall.  Satellite  measurements  of 
thermal  radiation  at  8 to  12  ju  m are  routinely  used  to  estimate  cloud  altitude.  The 
measured  radiances  are  converted  to  the  temperatures  of  a blackbody,  weighted 
over  the  spectral  interval  of  the  sensor,  and  along  with  a temperature  altitude  rela- 
tion, used  to  estimate  cloud  altitude.  Colder  temperatures  correspond  to  higher 
clouds,  while  warmer  terr\peratures  correspond  to  lower  clouds.  The  altitudes  are 
often  underestimated  when  satellites  view  cirriform  clouds  since  theory  and  ex- 
periment have  shown  that  cirriform  clouds  may  be  semitransparent  to  terrestrial 
radiation.  Reflected  sunlight  and  thermal  radiation  have  been  combined  to  infer 
cloud  types  by  Booth1®  and  by  Shenk,  Holub  and  Neff.  11  Clouds  which  are  highly 
reflective  and  warm  are  low  clouds  such  as  stratus.  Clouds  which  are  poorly 
reflective  and  cold  are  thin  high  clouds  such  as  cirrus.  Clouds  which  are  highly 
reflective  and  cold  are  thick  clouds  such  as  cumulonimbus  or  nimbostratus,  which 

are  expected  to  be  the  most  erosive  to  reentry  vehicle  nosetips.  Reflected  sun- 

12 

light  has  been  related  to  cloud  thickness  by  Reynolds  and  Vonder  Haar,  Griffith 
13  14 

and  Woodley,  and  Kaveney,  Feddes,  and  Liou  and  they  conclude  that  the  thicKes* 

clouds  appear  the  brightest.  Reflected  sunlight  and  also  thermal  emission  have 

been  related  to  rainfall  in  a number  of  investigations.  Martin  and  Scherer15  have 

1 6 

summarized  work  up  to  1973.  Griffith  et  al  have  tested  the  concepts  extensively 
using  visible  and  infrared  imagery  from  geosynchronous  satellites.  They  have 
found  that  areas  of  active  convection  and  rainfall  in  the  tropics  are  brighter  or 
colder  on  the  satellite  visible  or  infrared  photographs  than  inactive  regions. 

In  summary,  other  studies  have  shown  that  clouds  which  are  highly  reflective 
at  visible  frequencies  and  which  appear  rather  cold  at  infrared  window  frequencies 
tend  to  be  high,  thick,  and  produce  precipitation.  Such  clouds  are  also  likely  to  be 
erosive  to  reentry  vehicle  nosetips,  so  that  visible  and  infrared  satellite  data  are 
expected  to  be  useful  for  estimating  erosion  parameters. 

1.5  Approach 

The  Satellite  Correlation  study  used  an  empirical  approach  to  relate  satellite 
measurements  to  erosion  parameters.  Aircraft  flights  were  made  underneath 
satellites  in  a variety  of  cloud  conditions.  Non-linear  equations  were  developed  to 
relate  infrared  temperatures  and  visible  luminances  observed  by  satellite  to  erosion 
parameter  estimates  based  on  aircraft  measurements  of  cloud  and  precipitation 
particles.  Satellite  data,  aircraft  data,  and  a case  study  are  described  in  Section 2. 

Studies  were  first  directed  toward  techniques  to  estimate  the  altitudes  and  ice 
content  of  cirrus  clouds.  The  studies  were  justified  by  early  observations  that 
cirrus  clouds  eroded  the  nosetips  of  high  velocity  reentry  vehicles  and  by  the  fact 

(Because  of  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 
See  Reference  Page  79,  for  References  10  through  16. ) 
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that  existing  climatologies  of  cirrus  cloud  were  quite  limited.  Aircraft  and  satel- 
lite measurements  taken  in  1974  demonstrated  that  cirrus  altitude  and  ice  content 
could  be  estimated  to  useful  accuracies  from  multispectral  infrared  radiances 
measured  by  the  NOAA  Vertical  Temperature  Profile  Radiometer  (VTPR).  How- 
ever, by  1975,  customer  emphasis  shifted  from  cirrus  clouds  to  all  clouds  with 
particles  larger  than  about  50  pm.  For  this  more  general  problem,  a combination 
of  visible  and  infrared  data  from  NOAA  Scanning  Radiometers  (SR)  were  found  to 

be  more  useful  than  VTPR  data,  because  the  VTPR  technique  was  saturated  by  ice 

_2 

contents  in  excess  of  about  100  gm  integrated  through  a vertical  column.  Pinal 

1 7 

VTPR  studies  were  reported  by  Hunting. 

Preliminary  correlations  were  developed  between  satellite  visible  and  111  data 
and  weather  erosion  as  represented  bv  the  Fnvironmental  Severity  Index  (ESI),  de- 
fined in  Section  2.  5.  The  correlations  were  applied  to  a sample  month  of  data  over 
Eurasia  ami  maps  showing  areas  exceeding  certain  threshold  values  were  generated. 
A time  lapse  film  showing  daily  distribution  of  ESI  values  greater  or  equal  to  two 
for  one  month  over  the  area  was  also  produced.  A successful  demonstration  of  this 
work  in  April  1976  was  followed  by  a request  for  a full  year  of  climatology  over 
Eurasia,  the  Kwa.jalein  and  Wallops  areas. 

Since  197  6,  the  preliminary  correlation  between  ESI  and  satellite  data  was 
supplemented  by  more  cases  of  satellite  and  aircraft  measurements  in  heavy 
weather  situations.  Some  cases  of  simultaneous  satellite  and  radar  measurements 
of  clouds  were  also  taken.  Correlations  were  developed  for  erosion  parameters 
other  than  the  ESI.  Calibrations  were  acquired  for  NOAA  geostationary  satellites 
so  that  the  correlations  based  on  polar  orbiting  satellites  could  be  applied  to  geo- 
stationary satellites,  which  provide  data  over  most  of  the  Western  Hemisphere  at 
30-min  intervals.  The  calibrated  geostationary  satellite  data  were  used  in  real- 
time to  forecast  conditions  for  erosion  measurements  at  missile  test  ranges. 

Further  details  of  the  observational  programs,  correlation  studies,  clima- 
tologies, real-time  applications,  verifications,  and  important  computer  programs 
follow. 

2.  OBSERVATIONS 
2.1  Salrlliti' Data 

Data  from  the  NOAA  1TOS  series  of  satellites  were  used  to  correlate  erosion 
parameters  and  satellite  observations.  These  correlations  were  then  applied  to  a 


17.  Hunting.  J.  T.  (1976)  Cloud  Properties  from  Satellite  Infrared  and  Visible 
Measurements.  AKcTl., -TR -7 6 -02&7 , 6 pp. 
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year  of  NOAA  archive  to  produce  climatologies  of  certain  erosion  parameters. 

The  NOAA  Scanning  Radiometer  (SR)  archive  was  the  most  suitable  data  base  for 
several  reasons.  First,  data  were  available  every  day  starting  in  April  1973. 
Therefore,  data  were  available  for  the  last  10  months  of  the  period  February  1973 
to  January  1974,  which  had  been  selected  by  ABRES  for  intensive  studies  of  weather 
over  target  areas.  Second,  the  data  were  available  for  the  entire  Northern  Hemis- 
phere, which  includes  areas  where  calibration  flights  were  made,  test  ranges,  and 
target  areas.  Finally,  the  data  were  mapped,  calibrated,  and  available  on  mag- 
netic tape. 

Data  from  the  geosynchronous  satellites  GOES  1 located  at  75°  W and  SMS  2 
located  at  135°  W were  used  in  real  time  to  plan  aircraft  sounding  flights  in  con- 
junction with  passes  of  the  polar -orbiting  NOAA  satellites,  to  produce  real  time 
aerial  distributions  of  erosion  parameters  in  the  Wallops  and  Kwajalein  test  ranges, 
and  to  verify  some  of  the  established  relationships. 

2.  1.  1 ITOS  SATELLITE  DATA 

Data  which  were  correlated  to  aircraft  measurements  were  from  Scanning 
Radiometers  (SR)  aboard  the  polar  orbiting  NOAA  ITOS  2,  3,  4,  and  5 satellites. 

The  one  year  climatology  was  from  NOAA  2 (April  1973  to  19  March  1974)  and 
from  NOAA  3 (20  March  to  31  March  1974). 

18 

Details  of  the  NOAA  SR  archive  have  been  described  by  Conlan.  The  archive 
includes  simultaneous  IR  and  visible  measurements  once  per  day  over  the  Northern 
Hemisphere.  The  satellites  are  sun-synchronous  and  pass  over  mid-latitudes  at 
about  0900  to  1000  local  time.  The  visible  sensors  respond  to  light  at  wavelengths 
0.52  to  0.72  pm  except  for  NOAA  5 whose  range  is  about  0.44  to  0.92  pm.  The  IR 
sensor  responds  at  10.4  to  12.  5 pm  in  the  IR  window  region.  Energy  calibrations 
are  applied  to  provide  a data  bank  as  homogeneous  as  possible.  Variations  in  the 
spectral  response  of  different  visible  or  IR  sensors  are  considered  Insignificant 
for  our  application.  Both  the  visible  and  IR  are  archived  in  an  8 -bit  format.  The 
visible  code  is  proportional  to  luminance  in  foot -Lamberts.  The  maximum,  or 
saturation  value  of  254  is  approximately  the  luminance  of  a perfectly  diffuse  or 
Lambertian  reflector  at  the  top  of  the  Earth's  atmosphere  when  the  sun  is  overhead. 
Throughout  this  study  the  visible  data  are  retained  in  the  8-bit  scale  of  0 to  254. 

The  coded  IR  values  are  converted  to  "K  according  to  Table  1 for  subsequent  use. 
The  archive  of  visible  has  been  normalized  for  solar  elevation  by  dividing  each  value 
by  cosine  9 where  9 is  the  zenith  angle  of  the  sun.  Observations  in  the  IR  array  are 
much  less  sensitive  to  viewing  geometry.  The  IR  archive  includes  a small  correc- 
tion for  limb  darkening  due  to  water  vapor  absorption.  A daily  quality  control  is 
used  to  monitor  the  visible  and  IR. 

18.  Conlan,  E.  F.  (1973)  Operational  Products  from  ITOS  Scanning  Radiometer  Data, 
NOAA  Technical  Memorandum  NEtlij  b'2,  57  pp. 
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Table  1.  Conversion  of  [R  Coded  Values  to  °K 
for  the  NOAA  Scanning  Radiometer  Archive 


Coded  Value 

Temperature  'K 

0 

164.  0 

1 

165.0 

78 

242.  0 

79 

242.  5 

80 

243. 0 

81 

243.  5 

254 

330.  0 

255 

Missing 

Data  from  12  or  13  consecutive  orbits  are  archived  in  the  form  of  2048X2048 

* 

arrays  for  each  hemisphere.  The  array  of  visible  and  IR  data  is  aligned  with  the 
conventional  numerical  weather  prediction  grid,  so  that  the  arrays  are  equally 
spaced  on  a polar  stereographic  projection  as  shown  in  Figure  1.  The  spacing  of 
adjacent  grid  points  on  the  surface  of  the  Earth  increases  from  about  6 km  at  the 
Equator  to  about  13  km  at  the  poles  as  shown  in  Figure  2.  The  grid  spacing  is 
1/32  of  the  spacing  of  the  National  Meteorological  Center  grid.  Horizontal  reso- 
lution of  the  original  satellite  data  is  roughly  3 km  for  the  visible  channel  and  about 
6 km  for  the  IR  channel  near  the  satellite  subpoint.  The  grid  is  in  some  instances 
larger  and  in  other  instances  smaller  than  the  horizontal  resolution  of  the  original 
radiometer  data,  depending  on  whether  the  data  are  visible  or  IR,  on  the  zenith 
angle  of  observation  by  the  satellite  and  the  latitude  of  the  grid  point.  When  the 
original  radiometer  data  are  of  finer  resolution  than  the  array,  for  example,  when 
the  radiometers  are  looking  straight  down,  no  attempt  is  made  to  average  values 
close  to  the  array  point— only  the  last  value  processed  is  retained. 


* 

This  archive  was  changed  to  1024  X 1024  on  15  September  1976.  At  times  of  addi- 
tional cloud  truth  observations  after  this  date,  satellite  passes  were  specially  pro- 
cessed for  us  on  the  original  scale  of  2048  X 2048  thus  assuring  the  best  possible 
data  for  correlation  or  verification  purposes. 

16 
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Figure  1.  Satellite  Archive  Data 
Array  Over  the  Northern  Hemisphere 


Figure  2.  Satellite  Archive 
Data  Mesh  Size  vs  Latitude 


In  spite  of  the  normalization  of  the  visible  for  the  solar  zenith  angle,  incon- 
sistencies in  brightness  are  often  noted  at  the  swath  edges  between  satellite  passes. 
These  inconsistencies  can  be  seen  in  Figure  3,  an  image  generated  from  archived 
data.  Where  the  satellite  looks  eastward  toward  the  morning  sun,  forward  scat- 
tering causes  the  cloud  to  appear  excessively  bright.  Where  the  satellite  looks 
westward,  the  weaker  back  scattering  causes  the  same  cloud  to  appear  relatively 

dark.  These  irregularities  were  reduced  using  the  bi-directional  reflectance  model 

19 

of  Sikula  and  Vonder  Haar.  They  developed  correction  factors  from  available 
summaries  of  cloud  reflectivities  observed  by  aircraft  or  satellites  at  different 


19.  Sikula,  G.  J.  , and  Vonder  Haar,  T.  H.  (1972)  Very  Short  Kange  Local  Area 

Weather  Forecasting  Using  Measurements  from  Geosynchronous  Meteorologi- 
cal Satellites,  l-'inalHeport  for  Contract  F 19628-7  1 -C- 007  3.  AFC’HL-72 
07ETT 
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viewing  angles  am)  solar  elevations.  The  original  tables  were  expanded  to  include 

satellite  zenith  angle  intervals  of  (iO  to  70  and  70  to  80  for  use  with  GOFS  data.  The 

additional  factors  were  based  on  data  summarized  bv  liaschke.  Yonder  Haar, 

20 

Pasternak,  and  Bandeen.  The  geometry  of  reflectance  is  defined  by  Figure  4 
with 

C the  solar  zenith  angle, 
il  I the  solar  azimuth  angle, 

(1  the  satellite  zenith  angle, 
l the  satellite  azimuth  angle, 

t p equal  to  the  relative  azimuth  angle. 

From  these  angles  correction  factors  can  be  found  in  Table  2.  An  observed  value 
of  visible  luminance,  or  brightness,  already  corrected  for  the  sun's  zenith  angle, 
is  multiplied  by  the  appropriate  correction  factor. 


Figure  3.  Visible  Sensor  Imagery  Generated  From  the 
Digital  Archive  for  14  and  lf>  February  1875.  Coast- 
lines do  not  appear  in  the  archive 

liaschke,  F.  , Yonder  llaar,  T .If.,  F 
The  Radiation  Balance  of  the  Forth 
lla.liation  Aloasuroments.  NASA  T> 


Atmosphere  System  from  Nin.nus  3 


satellite 


Table  2.  Correction  Factors  (x  / r)  to  Adjust  Observed  Brightness  to  Reflectance 
from  an  Isotropic  Surface  With  the  Sun  Overhead 


2-A 

2-B 

2-C 

0‘S(<  25”  \/r  = 1.  02  = constant 

25”  - ? < 45”  (Linear  interpolation  between  Tables  2-A  and 

45”  < C < 60° 

0” -30”  30” -60°  60° -90°  90”-120”  120”-150° 

2-C) 

150” -180 

O’  -10 

1.  00 

1.00 

1.00 

1.  00 

1.00 

1.  00 

10° -20° 

0.93 

0.  93 

0.  93 

0.  95 

0.  95 

0.  95 

20” -30° 

0.90 

0.  90 

0.93 

0.  93 

0.  93 

0.  90 

30” -40 

0.  87 

0.  87 

0.  90 

0.  93 

0.88 

0.  83 

40” -50” 

0.  83 

0.  83 

0.  85 

0.  85 

0.  80 

0.70 

50” -60” 

0.75 

0.75 

0.76 

0.76 

0.  68 

0.  53 

60” -70” 

0.73 

0.73 

0.  81 

0.77 

0.  62 

0.  50 

70”-8 0” 

0.  70 

0.  72 

0.  80 

0.75 

0.  56 

0.  39 

2-D 

60*  < C < 80” 

0”  -30” 

30” -60” 

60” -90” 

90” -120” 

120” -150” 

150” -180” 

0”-10” 

1.  04 

1.  04 

1.  02 

1.  00 

1.  00 

1.  00 

10” -20” 

1.  07 

1.04 

1.  00 

1.  00 

0.  98 

0.  98 

20” -30” 

1.  04 

1.04 

1.00 

0.  97 

0.93 

0.  93 

30” -40” 

0.  97 

0.97 

0.97 

0.  93 

0.  86 

0.  80 

40” -50” 

0.  84 

0.  87 

0.  92 

0.  86 

0.74 

0.70 

50” -60” 

0.72 

0.78 

0.  82 

0.75 

0.  64 

0.  60 

60” -70” 

0.  68 

0.78 

0.83 

0.74 

0.50 

0.  41 

e 

O 

CO 

1 

e 

O 

0.  62 

0.74 

0.79 

0.  68 

0.  39 

0.  28 

19 


In  order  to  apply  correction  factors,  the  three  bi-directional  angles  must  be 
calculated.  It  was  decided  to  do  this  for  each  10  X 10  data  array  in  the  archive. 

This  corresponds  to  a square  of  about  100  km  on  a side  at  temperate  latitudes.  This 
increment  of  resolution  provides  angles  within  a degree  or  two  of  being  correct  for 
all  points  in  the  SR  array.  Inspection  of  Table  2 shows  that  angular  errors  of  this 
magnitude  produce  correction  factor  changes  that  are  considered  well  within  the 
overall  accuracy  of  the  table. 

The  hemisphere  archives  do  not  list  which  array  points  belong  to  a particular 
orbit;  therefore,  it  is  necessary  to  make  this  determination  first.  Then  the  three 
angles  for  the  center  of  each  10  X 10  array  can  be  calculated  knowing  the  ephemeris 
data  of  the  satellite  and  sun.  An  example  of  array  determination  by  orbit  for  the 
Northern  Hemisphere  for  14  and  15  February  1975  is  given  in  Figure  5.  This  pat- 
tern identifies  the  orbits  of  the  composite  image  shown  in  Figure  3. 

In  Figure  5,  orbit  1134  represents  the  earliest  data  in  the  array,  while  orbit 
1146  represents  the  latest  data,  which  is  12  orbits  or  23  hr  later.  The  pattern  of 
orbits  does  not  exactly  exhibit  circular  symmetry,  due  to  the  method  of  entering  the 
orbital  data  on  the  array.  All  of  the  data  for  a given  orbit  is  written  on  the  array 
when  the  orbit  is  processed;  however,  data  from  orbits  processed  later  will  re- 
place the  data  processed  earlier  when  a spatial  overlap  of  orbits  exists.  When  all 
of  the  data  for  a particular  orbit  is  entered  in  the  archive,  the  pattern  looks  like 
orbit  114  6,  which  was  in  fact  the  last  orbit  processed  so  that  its  eastern  side  re- 
placed some  of  the  data  for  orbit  1145  while  its  western  side  replaced  some  of  the 
data  for  orbit  1134.  When  orbits  are  processed  in  order,  the  pattern  looks  like 
orbits  1137  to  1145  on  Figure  5.  When  orbits  are  processed  out  of  order,  the 
pattern  looks  like  orbit  1135,  which  was  written  on  the  array  at  some  later  time 
than  orbit  1136,  and  replaced  some  of  the  data  on  the  eastern  half  of  1136.  This 
reversal  of  orbit  processing  can  occur  when  two  orbits  of  data  are  stored  on  the 
satellite  tape  recorders  and  read  out  later. 

The  order  in  which  orbits  were  processed  was  determined  by  looking  at  the 
pattern  of  discontinuities  on  images  of  the  arrays.  In  general  the  pattern  remains 
similar  from  day  to  day  except  for  rotation  around  the  pole  and  occasional  varia- 
tions near  the  dateline,  at  about  308  E. 

The  size  of  the  satellite  data  array  which  was  selected  for  comparison  with 
cloud  truth  information  was  based  on  several  factors.  It  was  noted  in  the  intro- 
duction that  at  the  start  of  the  study  the  satellite  VTPR  data  were  of  principal  in- 
terest. These  data  are  for  square  areas  of  about  70  X 70  km  at  satellite  nadir 
enlarging  to  rectangular  areas  about  70  X 108  km  at  left  and  right  limits  of  the 
scan.  Through  most  of  the  data  gathering  period,  the  size  of  the  SR  sample  was 
made  roughly  the  same  as  the  size  of  the  VTPR  measurement  over  the  area  of  in- 
terest for  comparative  purposes.  This  size  area  also  proved  to  be  representative 
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of  the  area  sounded  by  aircraft.  Ideally,  cloud  truth  should  apply  to  vertical  col- 
umns throughout  the  cloud  and  be  instantaneously  determined  at  the  precise  time  of 
the  satellite  scan.  In  practice,  spiraling  ascents  or  descents  at  300 to 400  m min-1 
and  about  35  km  in  diameter  proved  to  be  the  fastest  flight  pattern  adequate  for 
particle  sampling. 


NOAA-4 


Figure  5.  Pattern  of  Satellite  Data  by  Orbits  for  the  Northern  Hemisphere 
14  and  15  February  197  5 


The  SR  data  are  routinely  placed  into  the  NMC  grid  with  care,  nevertheless 
small  position  errors  can  usually  be  found.  In  order  to  minimize  these  errors,  and 
place  the  data  as  closely  as  possible  over  the  aircraft  sounding,  printouts  of  the  SR 


area  of  interest  were  made.  Transparencies  showing  coast  and  latitude  — longitude 
tines  were  prepared  in  the  same  scale  as  the  printouts.  These  were  generated  from 
an  archive  of  coastlines,  rivers,  lakes,  and  geopolitical  boundaries  referenced  to 
latitude  — longitude.  * The  SR  printouts  were  then  analyzed  to  determine  the  shifts 
necessary  to  fit  the  landmarks.  In  cases  without  landmarks  near  the  area  of  inter- 
est, Very  High  Resolution  Radiometer  (VHRR)  images  obtained  from  NESS  were 
hand  gridded  by  anchoring  the  grid  to  visible  landmarks.  Positions  of  prominent 
cloud  features  were  determined  and  transferred  to  the  transparencies.  The  SR 
printout  was  then  shifted  to  the  best  fit.  Often  hoth  techniques  were  applied.  The 
average  grid  shift  to  the  SR  data  over  locations  in  the  eastern  US  was  up  1,  left  1 
data  point.  The  maximum  shift  found  was  three  data  points. 

Neither  of  these  techniques  could  be  used  for  cases  near  Kwajalein  in  the 
Pacific  Ocean  because  of  the  lack  of  landmarks  and  VHRR  imagery.  However, 
through  the  use  of  photographs  at  )6-min  intervals  of  the  WSR-57  radar  PPI  scope, 
small  isolated  echoes  of  showers  could  be  pin-pointed.  Their  corresponding  images 
could  often  be  pin-pointed  in  the  satellite  data  to  position  the  SR  array  accurately. 
This  approach  requires  the  radar  films  to  be  timed  correctly,  which  was  not  always 
the  case.  Furthermore,  by  checking  the  location  of  towers  and  the  atoll  outline 
during  a period  of  anomalous  propagation,  it  was  found  that  the  azimuth  scale  on 
the  PP!  photography  required  a -2°  correction. 

IHie  to  the  variable  nature  of  tropical  clouds  and  showers,  sample  SR  areas  had 
to  be  reduced  to  a few  data  points  in  some  cases.  In  the  extreme  case,  the  satellite 
sample  had  to  be  reduced  from  the  standard  70  X70  km  (about  50  data  points  at 
35”  N to  as  little  as  7. 5 X 7.  5 km  (one  data  point  at  9'  N).  None  of  the  Kwajalein 
cases  were  included  in  the  correlations,  but  they  are  discussed  in  Section  5 on 
Verification. 

Satellite  variables  which  were  extracted  from  the  area  through  which  the  air- 


craft sampled  were  as  follows: 

(1)  ITT.  average  IR  temperature,  *K, 

(21  Hjpp  altitude  in  km  corresponding  to  TTT, 

(31  ci  TTT.  standard  deviation  of  the  IR  temperature. 
(41  o IR  (ITT,  IR  coefficient  of  variation, 


(51  TT^,  average  normalized  visible  brightness  in  arbitrary  units, 

(HI  a B,  standard  deviation  of  brightness  normalized  only  for  sun  zenith  angle, 
(71  cr  B/-r.  brightness  coefficient  of  variation. 

*The  original  archive  was  prepared  by  NESS;  it  was  refined  to  one  point  every 
one -tenth  of  a degree  for  some  areas  by  Cozzens.-l  The  refined  archive  is  now 
used  by  NESS  for  all  gridded  displays. 

21.  Cozzens.  D.  F.  (19741  Technique  to  Earth  I.ocate  and  Analyze  Oata  from 

Satellite  Sensors,  pp  13-3(5.  Burke,  L.  (19741  Met  and  ^laih.  Anal.  Pinal 
Rept.  Regis  Uollcge,  AFCRL-TR -74 -0558, 
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The  above  variables  were  entered  in  a stepwise  multiple  repression  propram 
to  determine  the  best  coefficients  to  predict  erosion  parameters  as  computed  by 
simultaneous  aircraft  soundings.  The  results  are  discussed  in  Section  3 of  this 
report. 

2.  1.2  OOES-I  and  SMS-2  GEOSYNCHRONOUS  SATELLITE  DATA 

Geosynchronous  satellite  data  were  not  available  to  the  Satellite  Correlation 
study  during  the  first  series  of  flights  in  the  winter  of  1973  and  1974.  Subsequently, 
our  first  use  of  the  geosynchronous  visible  and  IR  imagery  was  for  planning  air- 
craft sounding  flights.  Sectorized  pictures  were  received  by  C-5  telephone  line 
from  NESS.  In  the  beginning  these  had  to  be  gridded  using  our  own  library  of  grids 
along  with  NASA  ephemeris  predictions.  Later  the  images  were  received  already 
gridded.  Their  use  helped  greatly  in  forecasting  sounding  sites  4 to  6 hours  in  ad- 
vance of  aircraft  flights.  A more  significant  use  of  geosynchronous  data  was  pos- 
sible after  the  Man-computer  Interactive  Data  Access  System  (McIDAS)  became 
operational  at  AFGL  since  digital  data  could  be  obtained  and  processed.  Once  pro- 
grams were  perfected  to  navigate  the  data  and  convert  it  to  NOAA  satellite  equiva- 
lents, it  became  possible  to  print  out  visible  and  IR  data  and  erosion  estimates 
over  small  specified  areas  in  real-time. 

General  information  on  the  geosynchronous  satellites  has  been  reported  bv 
■>  •> 

Bristor““  and  is  summarized  here.  Visible  sensors  respond  in  the  0.55  to 
0.75  pm  range  while  IR  sensors  respond  in  the  10.5  to  12.5  pm  range.  Their 
spectral  response  characteristics  are  reasonably  close  to  those  of  the  NOAA  series 
of  satellites  thus  facilitating  the  conversion  process.  Before  reaching  McIDAS  the 
geosynchronous  data  are  "stretched"  so  that  the  rate  of  data  transmission  is  slower 
and  adjustments  are  applied  to  the  visible  to  make  each  scan  line  compatible  with 
adjacent  lities.  The  data  stream  is  then  returned  to  its  respective  satellite  and  re- 
transmitted. Information  from  GOKS-1,  located  at  0*,  75"  W,  and  SMS -2  over  0°. 
135“  W can  be  received,  but  not  simultaneously,  at  AFGL. 

In  order  to  compare  data  from  the  polar-orbiting  NOAA  series  of  satellites  with 
data  from  the  geostationary  satellites  GOES-1  and  SMS-2,  available  information  on 
sensor  calibration  was  reviewed.  No  conversion  factors  were  necessary  for  IR  data 
since  all  the  satellites  had  IR  sensors  calibrated  in  “K.  For  visible  data,  ground 
calibrations  established  before  lautu  h were  used  to  develop  the  following  equations. 

To  convert  8-bit  visible  data  from  the  geostationary  satellites  to  the  8 -bit 
NOAA  satellite  scale: 


Bristor,  C.  L.  (19751  Central  Processing  and  Analysis  of  Geostationary  Satellite 


Data,  NOAA  Tech  Mem. 


Washington. 
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(1) 

(2) 


where  CL,  is  the  8 -bit  NOAA  count, 

N 

r . is  the  6 -bit  GOES-1  count,  and 
G 1 

CS2  is  the  6 -bit  SMS -2  count. 

Brightness  normalization  for  solar  elevation  is  performed  by  McIDAS  as  follows 


B 


C 


N 


N ” cos  9 ' 


(3) 


where 

cos  9 = sin  6^,  X sin  6g  + cos  6^,  X cos  6g  X cos  (A^,  - Ag)  , (4) 

with  the  symbols  6 and  A referring  to  latitude  and  longitude,  and  the  subscripts  C 
and  S referring  to  cloud  and  sun. 

Corrections  for  bi-directional  reflectance  (Table  2)  cannot  be  calculated  pre- 
sently within  McIDAS  due  to  limited  storage  capacity.  These  are  currently  calcu- 
lated externally  and  then  a single  factor  is  inserted  in  McIDAS  to  produce  fully 
normalized  brightness  for  the  area  of  interest. 

Maximum  resolution  at  the  sub -point  of  the  geostationary  satellites  is  8 km  in 
the  IR  and  1 km  in  the  visible;  therefore,  the  visible  should  be  combined  to  provide 
a value  comparable  to  the  smallest  possible  IR  element.  Within  the  8 km  IR  element 
there  are  eight  lines  of  visible  each  containing  eight  pixels.  Present  practice  is  to 
average  all  eight  pixels  out  of  one  of  the  eight  lines  and  to  assign  the  value  to  the 
entire  square.  All  eight  lines  cannot  be  averaged  due  to  computer  limitations. 

Actual  sizes  of  one  of  these  IR  elements  are,  for  example,  over  Wallops  Island 
9.  8 km  N-S  X 6.  2 km  F-W  (from  GOES-1);  and  over  Kwajalein  8.  1 km  N-S  X 15.  2 
km  E-W  (from  SMS-21. 

The  navigational  accuracy  of  positions  at  sea  level  depends  on  how  recently  an 
orbital  perturbation  was  introduced  and  how  close  the  navigational  landmarks  are 
to  the  area  of  interest.  For  example  when  a rather  distant  landmark  like  Lake 
Titicaca  in  S.  America  is  used  to  navigate  in  the  central  USA,  an  accuracy  of 
2 to  4 km  can  be  attained  after  four  pictures.  After  8 to  10  pictures  the  accuracy 
improves  to  about  1 km. 
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It  is  important  to  remember  that  when  zenith  angles  to  the  satellite  *re  • i 
to  large  the  position  of  high  clouds  is  in  considerable  error.  For  exar  pic. 
Kwajalein,  where  the  zenith  angle  to  SMS-2  is  65°,  a cloud  10  km  high  * t • • 

gated  about  20  km  west  of  its  true  position  unless  corrections  are  applie  i. 

2.2  Aircraft  Data 

Aircraft  flights  in  clouds  viewed  bv  satellites  were  made  by  four  liffer  «•>  • > ■ 
craft  from  January  1074  toOctober  1076.  A USAF  C130-A  aircraft  read* 
soundings,  while  a C130-E  made  15,  a Cessna  Citation  operated  bv  MR  1, 
two  and  a l.ear  36  operated  by  Aeromet,  Inc.  and  MRI,  Inc.  made  three, 
aircraft  and  Citation  were  used  until  2 March  1076  and  the  l.ear  let  w.i- 
that  date. 

Aircraft  flights  were  directed  to  suitable  cloud  conditions  by  i lou  l • ■ 
from  AFGL.  The  forecast  was  primarily  a 2-  to  6-hr  extrapolation  of  i».'l 
pictures,  hourly  w-eather  station  reports  and  NWS  facsimile  products  t ••  • 
satellite  pass.  The  forecasts  were  made  necessary  by  time  required  to  • 

airspace  from  the  FAA  and  to  fly  to  the  area  of  interest.  Early  in  the  t>r 
flights  were  locked  into  the  forecasted  location.  Later,  however,  the  air  • • 
acquired  some  ability  to  make  last-minute  changes  in  the  flight  deslu 
on  their  own  observations. 

2.2.  1 C-130  AIRCRAFT  AND  DATA  REDUCTION 

The  C-130  aircraft  were  based  at  Langley  AFB  near  Norfolk.  Virgo  < 
ings  were  made  as  far  as  1000  km  inland  from  Langley  but  not  over  sever  ■ 
kilometers  out  to  sea.  Most  sounding  patterns  were  spiral  descents  ar  ■ . 

gational  beacons  with  entry  into  the  highest  clouds  timed  to  coincide  c-  be-- 
with  the  satellite  pass.  Since  airspace  Suitable  for  rapid  ascents  ami  ie~  . 
difficult  to  find  over  the  eastern  U.  S.  , it  w as  often  advantageous  to  plan 
zones  restricted  to  military  flights.  The  C-130  aircraft  ceiling  is  about  ! 
a few  cases  cloud  tops  exceeded  this  level  so  top  altitudes  had  to  be  estr  at« 
the  flight  director. 

Additional  eye  observations  from  the  aircraft  were; 

(a)  The  depth  of  cloud  required  to  obscure  the  sun's  disk. 

(b)  Halo  phenomena, 

(cl  Ice  particle  form,  size,  and  intensity  as  seen  on  the 
"snowstick"  probe, 

(dl  Presence  of  rime  and  its  rate  of  formation, 

(e)  Change  of  particle  phase, 

(f)  Cloud  levels,  types,  fractional  coverage,  and 

(g)  Estimates  of  cloud  optical  density. 
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Hydrometeor  and  atmospheric  probes  aboard  the  C- ISO's  are  listed  in  Figure G 
along  with  their  depiction  on  the  aircraft.  Unfortunately  data  from  the  hydro- 
meteor probes  frequently  had  limited  utility  due  to  malfunction,  lack  of  a full  under- 
standing of  what  the  instrument  sensed  or  difficulties  in  performing  an  analysis. 


1 OFW  COIN T HYGROME  TEH  PROBE 

2 PMS  1-0  AXIAL  SCATTER  PR0B£(2-3O|i) 

3 RIS  1-0  PRECIP.  PROBE  (300- 4 900 u ) 
4.HYDR0METE0R  TOIL  SAMPLER 

9 PMS  I D CLOUO  PROBE  (20-300p) 

6 TOTAL  AIR  TEMPERATURE  PROBE 
T EWER  PROBE 

B PMS  2-0  PRECIP  PROBE (200- 6400 p) 
9 PMS  2-0  CLOUD  PROBE  (29-800(1) 
■OJW  CLOUO  WATER  probe 
II  POP  8/E  COMPUTER  B LINE  PRINTER 


12  EORMVAR  HYOROMETEOR  REPLICATOR 
IS  VISUAL  HYOROMETEOR  PROBE 


MINS  9 DOPPLER  RADAR 
13  AN/APOI22  Kg  S 9CM  WEATHER  RADAR 
16  16mm  NOSE  CAMERA 


1 7 PROBE  LIGHT 


Figure  G.  Atmospheric  Probes  Aboard  t he  C-130K  Aircraft 
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Analysis  of  the  hydrometeor  foil  samples  (Church,  Poos,  and  Spatolal  ' and  the 
formvar  replicator  were  particularly  difficult.  On  the  other  hand,  l.WC  profiles  of 
useful  accuracy  could  he  constructed  when  the  available  data  were  coupled  with  the 
ex<  ellent  eye  observations,  surface  observations  of  precipitation  rate,  and  PPI 
radar  data. 

When  precipitation  reached  the  ground,  a rainfall  rate  over  the  sample  area 
was  estimated  from  rain  gauge  and  radar  data.  The  precipitation  rate  was  con- 
verted to  mass  density  using  the  logarithmic  relationships  for  rain,  small  snow, 

and  large  snow  shown  in  Figure  7.  These  relationships  were  based  on  summaries 

24 

by  Plank.  The  estimated  mass  density  for  precipitation  was  simply  extrapolated 
upward  to  the  generating  altitudes  of  the  clouds.  The  mass  density  estimated  for 
rain  was  increased  to  a consistent  value  for  snow  above  the  melting  level. 


Figure  7.  Precipitation  Rate  — Mass 
Density  Relationships  for  Rain, 
Small  Snow,  and  Large  Snow 


0 1 I 10  100 

PRECIPITATION  HATE,  mm  hr  ' 

For  each  C-130  sounding,  an  altitude  profile  of  hydrometeor  density  was  con- 
structed which  was  as  consistent  as  possible  with  all  the  sources  of  data.  Recent 
measurements  of  ice  mass  vs  temperature  bv  the  University  of  Washington  and 


23.  Church,  .1.  F.  , LtOol,  Poes,  K.K.,  and  Spatola,  A.  A.  ( lit? SI  The  Continuous 

Aluminum-Foil  Hydrometeor  Sampler;  Design,  Operation,  Data  Analysis 
Procedures,  and  Operating  Instructions,  ArCrtl  -"tR -?Ti -M7n.  7 0 pp. 

24.  Plank,  V.  (5.  (13741  Hvdromctoor  Parameters  Determined  from  the  Radar  Data 

of  t lie  SAMS  Rain  Frosion  f’rogram,  SAMS  flepi.  Vo.  77  A JY'R I -TR  -?  1 - 

o24<i.  flft  pp. 


27 


25 

Meteorology  Research  Inc.  which  are  given  by  Hobbs,  Racike  and  Atkinson  and 

2 6 

profiles  through  16  storms  at  Wallops  Island  as  shown  by  Berthel,  were  used  as 
basic  guidance  in  constructing  profiles  of  liquid  and  ice  mass.  The  expertise  of 
Dr.  Robert  M.  Cunningham,  Senior  Scientist  at  AFGL,  was  drawn  upon  to  finalize 
all  profiles.  An  example  of  a sounding  and  the  construction  of  a profile  is  given 
in  Section  2.  3. 

2.2.2  LEAR  36  ,1ET  AIRCRAFT 

The  Lear  jet  ceiling  is  about  13.7  km;  all  flights  in  the  USA  topped  the  highest 
clouds.  Its  sounding  pattern  was  a racetrack  with  a length  of  about  50  km  along  a 
radial  from  a navigational  aid.  Climbing  soundings  were  usually  made  to  prevent 
fogging  of  the  PMS  probes. 

Details  of  the  High  Altitude  Weather  Aircraft  Data  System  (HAWADS)  are  given 
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bv  Hooker  and  Windes.  Instrumentation  is  essentially  the  same  as  that  aboard 
the  C-130's  except  for  the  lack  of  the  formvar  replicator,  foil  sampler,  Doppler 
radar,  and  EWER  probe.  In  most  cases  the  two-dimensional  (2D)  PMS  probe  could 
be  used  to  classify  particle  type  as  input  to  the  mass  computations  made  from  the 
one -dimensional  (ID)  data. 

The  greatest  remaining  difficulty  was  in  making  allowance  for  the  horizontal 
motion  of  the  hydrometeor  system  during  the  aircraft  sounding  in  order  to  develop 
a LWC  profile  at  the  exact  time  of  the  satellite  scans.  Radar  precipitation  patterns 
showing  echo  velocities  and  trends  were  the  principal  source  of  information  used 
to  make  these  adjustments. 

In  general  the  Lear  system  benefitted  from  earlier  experiences  with  the  C-130 
systems  and  the  Lear  soundings  were  an  improvement  over  the  C-130  flights  be- 
cause of  improved  instrumentation  and  understanding  of  the  measurements,  greater 
flexibility  in  site  location  and  less  time  consumed  en  route  to  the  site. 

2.3  Kiular  Data 

2.3.1  SPANDAR  RADAR 

The  SPANDAR  (Space  Range  Radar)  is  located  at  Wallops  Island,  Virginia.  It 
is  a high  resolution  (0.48’  beam  width)  "S"  band  (10.7  cm  wavelength)  tracking 
radar  having  a peak  power  output  of  3 MW.  Minimum  detectable  mass  for  bullet - 
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rosette  particles  is  about  0.001  g m at  a distance  of  50  km.  The  radar  is  cali- 
brated by  measuring  echo  returns  directly  over  a disdrometer. 

A program  was  initiated  to  gather  data  from  this  radar  at  times  of  NOAA  satel- 
lite passes  and  later  at  times  of  GOES  1 satellite  scans  over  the  area.  Radar  scans 
were  taken  over  a 34°  azimuth  sector.  These  data  were  to  be  processed  to  give 
integrated  mass  over  squares  corresponding  to  NOAA  grid  squares.  The  program 
for  analysis  of  four  cases,  three  ot  which  contained  detectable  echoes,  was  not 
completed  in  time  for  this  report. 

Data  from  SPANDAR  Range  Height  Indicator  (RHI)  scans  can  be  processed  by  a 
Liquid  Water  Content  Analyzer  to  estimate  a mass  profile  with  250  m increments 
of  height  and  also  the  ESI  . Radar  reflectivities  are  converted  to  mass  estimates 
by  means  of  logarithmic  relationships  based  on  radar  to  aircraft  data  comparisons 
for  various  particle-types  and  temperatures.  These  profiles  are  for  a column 
about  250  m along  the  beam  and  of  varying  width  across  the  beam.  At  a distance  of 
60  km  the  width  is  about  500  m.  Sample  volumes  are  therefore  very  small  com- 
pared to  samples  by  satellite,  nevertheless,  some  comparisons  of  ESIp  derived  by 
radar  and  satellite  are  discussed  in  Section  5,  Verification. 

2.3.2  TRADEX  RADAR 

The  TRADEX  (Target  Resolution  and  Discrimination  Experiment)  radar  is 

located  on  the  Kwajalein  Atoll  in  the  Pacific  Ocean.  It  is  also  an  "S"  band  radar 

having  a wavelength  of  10.  17  cm,  0.  3"  beam  width  and  minimum  detectable  mass 
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of  0.007  g m for  bullet -rosette  particles  at  50  km  distance. 

On  one  occasion  RHI  data  were  collected  from  this  radar  near  the  time  of  satel- 
lite scans,  and  a discussion  of  the  analysis  is  given  in  Section  5,  Verification. 

2.4  Sample  Cane— 8 October  1976 

Each  case  used  to  develop  the  general  relationships  between  satellite  and 
erosion  parameters  became  a small  study  in  itself.  All  cases  suffered  from  prob- 
lems; timing  differences  between  aircraft  and  satellite  observations  were  inevitable 
and  sometimes  there  were  missing  or  conflicting  measurements.  This  case  illus- 
trates good  aircraft  data,  some  timing  difficulties  and  an  unusual  problem  with  the 
satellite  data. 

2.4.1  PRE-FLIGHT 

During  late  afternoon  of  7 October  197  6 the  Lear  jet  aircraft  crew  was  alerted 
for  a possible  sounding  flight  from  Hanscom  AFB  at  Bedford,  Massachusetts  in  the 
New  York  state  or  New  England  area  on  the  following  day  at  1520Z  (10:20  AM  EST), 
the  time  of  the  NOAA -5  satellite  pass.  The  National  Environment  Satellite  Service 
in  Suitland,  Maryland,  was  notified  that  we  might  need  the  NOAA  pass  over  that  area. 


At  this  time  the  NOAA  archive  had  been  reduced  so  it  was  necessary  to  request 
special  processing  of  the  pass. 

On  the  morning  of  the  8th  it  was  evident  from  the  overnight  sequence  of  geo- 
synchronous satellite  pictures  that  the  thickest  clouds,  which  were  desired,  were 
moving  rapidly  northeastward  and  out  of  the  southern  sections.  It  was  therefore 
decided  to  go  to  Burlington,  Vermont  (BTV),  near  the  Canadian  border,  only 
30  min  distant  in  flying  time.  The  area  was  suitable  for  descent  to  low  levels  and 
hopefully  would  still  be  within  the  zone  of  thick  clouds  at  pass  time.  An  FAA  clear- 
ance was  obtained  for  an  ascent  and  descent  in  a racetrack  pattern  18  to  7 5 km  out 
of  BTV  on  the  180°  radial.  Plans  were  to  fly  from  Hanscom  AFB  to  Burlington, 
execute  a missed  approach  to  the  field,  then  climb  in  the  racetrack  pattern  to  cloud 
tops.  Breakout  would  be  slightly  ahead  of  pass  time  at  152 1Z.  Photographs  of  the 
top  would  be  made,  the  aircraft  would  reenter  the  tops  at  pass  time  and  descend  in 
the  racetrack  pattern  to  2 km,  minimum  safe  altitude  over  the  racetrack. 

NESS  was  notified  that  we  definitely  wanted  the  NOAA  pass  processed  in  the  old 
2048  X 2048  element  Northern  Hemisphere  format. 

2.  4.  2 GENERAL  SYNOPTIC  SITUATION 

Surface  fronts  and  isobars  are  shown  in  Figure  8 superimposed  on  the  digital 
visible  NOAA-5  satellite  data  which  was  gridded  and  converted  to  an  image  by  NESS 
through  the  use  of  an  arbitrary  gray  scale.  Also  indicated  on  Figure  8 are  the  posi- 
tions of  Burlington  (BTV),  the  weather  radar  station  at  Plattsburg  (PBG),  and  Albany 
(ALB),  the  nearest  RAOB  station.  The  area  of  the  racetrack  sounding  is  indicated 
by  arrows. 

In  the  sounding  area  cool  air  flowing  from  the  north  was  only  0.  5 km  deep;  a 
transitional  zone  between  cool  and  warm  air  was  found  between  0.  5 and  1.7  km. 
Above  this  level  winds  were  generally  from  the  SW,  averaging  72  km  hr  * at  the 
500  mb  level.  The  tropopause,  according  to  the  12Z  ALB  RAOB  was  at  12.7  km 
altitude.  The  PBGCPS-9  radar  showed  light  echoes  over  the  sounding  area  during 
ascent  and  before  pass  time,  but  no  echoes  at  or  after  pass  time.  Felloes  were 
moving  from  220"  at  88  km  hr  *.  Figure  9,  in  which  the  ground  pattern  has  been 
removed,  shows  echoes  north  and  west  of  the  racetrack  at  pass  time.  Aircraft 
sensors,  however,  detected  very  light  rain  until  termination  of  the  descent  at  2 km 
and  29  min  after  pass  time.  The  ceiling  at  BTV  was  400  to  500  ft  with  light  rain 
throughout  the  sounding  period. 


Figure  8.  Display  of  NOAA-5  Digital  Visible  Data  1512Z 
8 October  1976,  Surface  Synoptic  Pattern  at  1500Z  and 
Location  of  Reference  Points 


Figure  9.  PPI  Display  of  Weather 
Fchoes  at  Satellite  Pass  Time. 
Outer  range  circle  is  111  km.  Tin 
aircraft  sounding  took  place  along 
the  heavy  arrowed  line 


2.4.3  HYDROMETEOR  MASS  PROFILES 


Particle  sizes  and  their  numbers  are  measured  primarily  by  Particle  Measur- 
ing Systems  (PMS)  probes.  These  probes  are  optical  array  spectrometers  which 

were  developed  by  Knollenberg.  Recent  studies  of  their  response  to  ice  and  snow 
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particles  can  be  found  in  Knollenberg.  ’ The  one-dimensional  (1-D)  probes 

count  numbers  in  various  ranges  of  length.  To  reduce  these  data  to  mass  it  is 

necessary  to  classify  the  shapes  of  the  particles  counted  and  assign  pre-determined 

conversion  factors  for  each  shape.  Programs  for  converting  these  data  to  mass  are 
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described  by  Serbagi,  Rodenheiser,  and  Kaplan;  Belsky,  Kaplan,  and  Roden- 
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heiser;'  and  Belsky,  Francis,  Kaplan,  and  O'Neil. 

The  classification  is  greatly  facilitated  bv  data  from  the  PMS  two-dimensional 
(2-D)  probes  which  can  be  processed  to  display  particle  shapes  in  two  dimensions 
as  shown  in  Figure  10.*  In  this  display,  time  proceeds  from  left  to  right  along 
each  horizontal  strip.  Each  particle  that  enters  the  field  of  view  causes  a vertical 
line  whose  length  corresponds  to  1.  28  mm  particle  size.  To  the  right  of  the  line, 
the  particle  shape  is  displayed  in  the  same  scale.  The  time  corresponding  to  one 
sample  strip  varies  depending  on  the  size  and  number  density  of  particles.  Not 
all  particles  are  printed  out  in  this  display  due  to  the  huge  numbers  counted;  how- 
ever, ample  data  are  available  to  detect  changes  in  particle  character.  Figure  10 
shows  sample  strips  during  ascent  over  the  temperature  range,  + 5“C  to  -40°C. 
Raindrops  appear  somewhat  oval  but  should  appear  round.  The  elongation  is  caused 
by  a slightly  incorrect  air  speed  in  the  processing  program.  Drop  size  at  a temper- 
ature of  + 5®C  was  about  0.  5 mm.  At  0°C  clumps  of  crystals  were  samples;  one 
(circled)  was  about  2.5  mm  across.  As  the  temperature  decreases  particle  size 
also  decreases.  At  -ICC  and  -20I,C  dendritic  crystals  (circled)  about  0.5  mm  and 
0.4  mm  across  were  respectively  sampled.  From  these  printouts,  particles  en- 
countered during  the  sounding  were  classified  according  to  temperature,  altitude, 
and  time  as  shown  in  Table  3. 

From  the  1-D  count  data  and  eye  classification  of  2-D  data,  average  mass  w'as 
computed  for  30  sec  intervals  or  the  equivalent  of  about  200  m layers  in  the  atmos- 
phere. 

Mass  density  values  were  then  plotted  against  height  to  give  profiles  as  shown 
in  Figure  11.  These  are  for  particles  greater  than  about  50  pm  in  size,  generally 
referred  to  as  "precipitation"  particles.  Significant  observer  notes  and  temperatures 
are  also  shown.  Mass  profiles  of  liquid  cloud  particles  assumed  to  be  less  than 
50  pm  in  diameter  and  measured  by  the  .lohnson-Williams  sensor  are  plotted 
separately  as  shown  in  Figure  12. 


Methods  to  eliminate  the  need  for  1-D  data  by  2-D  counting,  sizing,  and  converting 
to  mass  estimates  were  under  study  at  the  time  of  this  flight  and  could  not  be  used. 


(Because  of  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 
See  Reference  Page  79,  for  References  28  through  32.  ) 
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Table  3.  Particle  Classification,  8 October  1076 


Particle  Classification  Change 

Temp.  (*C) 

Alt  (kin) 

Time(Z) 

Rain  to  large  snow 

+ 2 

3.  6 

143746 

I.arge  snow  to  small  snow 

-10 

5.  6 

144317 

Small  snow  to  bullet -rosettes 

-26 

7.9 

1447  16 

Bullet -rosettes  to  clear 

-42 

10.  0 

145100 

Clear  to  bullet -rosettes 

-45 

10.  2 

152701 

Bullet -rosettes  to  small  snow 

-29 

8.  4 

153317 

Small  snow  to  large  snow 

-24 

7.7 

153447 

I.arge  snow  to  rain 

+ 1 

3.  8 

154517 

Figure  11.  Vertical  Profiles  of  Mass  for  Particles  Greater  than  SO  |i  in  in  Diameter 
During  Ascent  and  Descent  8 October  107  6.  Temperatures  in  °C  are  shown  on  the 
left  with  observer  notes  on  the  right.  Points  are  PMS-1P  system  measurements. 
Points  with  crosses  (X)  are  averages  over  more  than  one  sampling  run 
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Figure  12.  Vertical  Profiles  of 
Mass  for  Particles  Less  than  50  pm 
in  Diameter  Measured  by  the 
Johnson-Williams  Sensor  During 
Ascent  and  Descent  8 October  107  6 


The  following  chronological  discussion  of  flight  observations  refers  mainly  to 
Figure  11  since  erosion  parameters  are  determined  from  the  profile  made  up  of 
precipitation  particles.  The  first  measurements  were  taken  in  a descent  into  BTV 
from  2.  1 to  0.  2 km  during  1421  to  1426Z,  almost  one  hour  before  pass  time.  The 
track  was  from  about  18  km  NW  of  BTV  to  BTV.  Light  to  moderate  rain  was  en- 
countered. Ascent  commenced  over  BTV  at  142630Z  over  a ceiling  of  about  125  m. 
The  aircraft  immediately  headed  to  the  racetrack  area  south  of  BTV  where  it  re- 
mained. During  ascent  the  lowest  clouds  were  topped  at  1.  8 km  as  noted  by  the 
observer  but  some  cloud  particles  were  measured  up  to  5.75  km  as  shown  by  Fig- 
ure 12.  The  observer  noted  another  liquid  cloud  layer  which  was  topped  at  3.  3 km. 
Very  light  rain  changed  to  very  light  snow  at  3.  6 km,  about  250  m below  the  freezing 
level.  A zone  of  agglomerates  averaging  1.0  mm  in  size  with  maxima  of  1.5  mm 
was  noted  at  5.  2 km.  This  agrees  well  with  a measured  maxima  of  mass.  Above 

8.  5 km  the  ice  cloud  was  not  solid  and  mass  decreased  significantly.  Mass  fell  to 
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less  than  1 mg  m v at  9.  65  km  and  at  145 1Z  the  aircraft  broke  through  wispy 

tops  at  10.0  km.  These  tops  were  therefore  2.7  km  below  the  tropopause  as  mea- 
sured near  ALB.  Character  of  the  cloud  tops  is  shown  in  Figure  13,  a photograph 
taken  at  9.  96  km.  The  aircraft  continued  to  climb  to  12.  96  km  then  it  descended 
to  reach  the  highest  clouds  at  1527Z,  a few  minutes  after  pass  time.  By  this  time 

the  ice  clouds  were  thinning;  patches  of  liquid  cloud  could  be  seen  through  thin  spots 
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as  shown  in  Figure  14.  Between  5.  35  and  4.  5 km  mass  dropped  below  1 mg  m in- 
dicating the  injection  of  drier  air  (from  the  SW)  and  gradual  weakening  of  the  cloud 
system.  A liquid  cloud  layer  was  entered  at  4.  5 km.  The  temperature  was  -2°C 
and  rime  built  up  rapidly.  Ice  crystals  were  found  and  they  grew  rapidly.  Snow 
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changed  to  rain  it  2.  8 km  atui  numerous  lavers  of  cloud  were  observed  down  to 
2.  1 km  where  the  sounding  terminated.  The  PMS  probes  fogged  up  at  the  2.  fi  km 
level. 


1 he  profile  at  pass  time  was  considered  to  be  the  average  of  t lie  ascent  and 
descent  profiles.  More  specifically,  average  values  were  determined  for  200  in 
lavers  centered  at  0.  1.  0.2  km.  etc,  , above  the  ground.  Two  sets  of  tabular  values 
are  made,  one  for  particles  above  about  50  ji  m in  si/e.  the  "precipitation  ' 'articles, 
and  one  which  includes  the  smaller  "cloud"  particles.  The  latte-  is  considered  to 
be  the  total  mass  or  1 WO, 


Figure  12.  Character  of  the  Cloud  Tops  Observed  from 
Top  l evel  of  10.0  km  During  Ascent 
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2.4.4  SATELLITE  DATA  FOR  8 OCTOBER  1976 

The  visible  satellite  data  were  received  and  processed  without  incident.  A 
visible  image  covering  the  area  of  interest  was  shown  earlier  in  Figure  8.  Digital 
data,  normalized  for  the  solar  zenith  angle,  are  shown  in  Figure  15.  The  aircraft 
racetrack  sounding  area  is  indicated.  It  was  assumed  that  the  cloud  pattern  re- 
mained unchanged  while  it  moved  at  the  velocity  of  the  radar  echoes  during  the 
sounding  interval.  Under  this  assumption  the  cloud  area  that  corresponded  to  the 
sounding  area  is  enclosed  by  the  parallelogram  shown  in  the  figure.  Ascent  took 
place  through  the  area  enclosed  above  and  to  the  right  of  the  aircraft  track  and 
descent  occurred  through  the  area  below  and  to  the  left  of  the  track.  The  IR  data 
from  this  pass  were  inadvertently  mapped  near  the  Aleutian  Islands.  Transforma- 
tion to  the  correct  area  involved  more  than  rotation  around  the  pole  because  of  scale 
changes  with  latitude  as  shown  earlier  in  Figure  2.  The  transformation  was 
achieved  by  printing  out  a large  area  of  visible  and  IR  between  two  points  of  positive 
identification.  The  IR  was  then  scaled  to  the  visible  and  tiie  correct  area  found. 
Coordinates  of  the  IR  are  therefore  different  from  the  visible  and  when  printed  out, 
as  in  Figure  16,  appear  rotated  approximately  90°  to  the  visible. 

* During  the  original  analysis  of  this  case  a different  echo  velocity  was  used  and 
therefore  a different  area  was  mapped  in  the  digital  data.  Since  the  satellite 
parameters  are  only  slightly  changed,  they  have  not  been  corrected  and  the 
regression  analysis  was  not  repeated. 
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Figure  15.  Aircraft  Sounding  Track  and  Area  of  Digital  Visible  Satellite 
Data  Assumed  to  Have  Passed  the  Track  During  Ascent  and  Descent. 
Areas  brighter  than  200  are  stippled 


The  contour  analysis  in  Figures  15  and  16  shows  that  relatively  cold  and  bright 
clouds  were  found  over  the  ascent  area  and  the  northern  end  of  the  racetrack.  This 
pattern  agrees  with  the  rain  pattern  as  seen  by  radar,  aircraft  observer  notes  and 
decrease  in  mass  from  ascent  to  descent. 
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Figure  16.  Same  as  Figure  15  Except  for  IR  Data  in  'K.  Areas  colder  than  -UP  ‘ 
stippled 


2.5  Summary  of  Satellite  ami  Erosion  Parameters  For 
All  Cases 

The  satellite  parameters  which  were  determined  for  each  case  were  les  t • 
in  Section  2.  1.  1.  Individual  values  for  each  case  are  given  in  columns  1 thr  u.  • 
of  Table  4.  Other  variables  listed  in  Table  4 were  derived  from  aircraft  ••  ei-  ,-  . 
ments.  They  are: 
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Column  No. 


i 

'< 


8 


;i 

10 


1 1 


12 


13 


I I 


AH 


I.WC 

i:si 

p 


The  total  thickness  in  km  of  all  cloud 

layers  having  liquid  or  ice  content  of 
, -3 

1 mg  m or  more. 

Total  liquid  and  ice  mass  for  the  column 
-2 

ingni  . 

Environmental  Severity  Index  for 
precipitation  particles  defined  as 


2 -3. 

ESIp  = J p 1 Id  1 1 (g  km  m ) 


(hi 


where  |i  is  the  density  of  hydrometeors 
(g  m ‘ ) and  11  is  the  altitude  above  the 
ground  (km). 

ESI„  Environmental  Severity  Index  for  all 

particles,  defined  in  the  same  way  as 

ESI  . 

P 

ADI  Accuracy  Degradation  Index. 

Ml  Range  error  in  feet. 

AS  Nose  tip  recession  in  inches. 


Values  of  ADI.  Ml.  and  \S  were  calculated  by  TRW,  Inc.  The  values  for  Ml 
and  AS  are  based  on  a TRW  erosion  model  for  a particular  design  of  reentry  system. 
The  calculations  use  the  altitude  profile  of  precipitation  mass  density  as  input. 
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I'.ililc  -1.  Dependent  Data  Sample1.  Satellite1  parameters  anil  e-oi  riaspeiiiilm^  mte^rateel 
elenul  mass,  elouei  thie’klie'ss,  anei  erosion  parameters 
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3.  RELATIONSHIPS  BETWEEN  SATELLITE  PARAMETERS  AND 
EROSION  PARAMETERS,  LIQUID  WATER  CONTENT 
AND  CLOUD  THICKNESS 


3.1  Transformation  of  Satellite*  Data  to  a Linear  Korin 

Relationships  between  the  satellite  data  and  simultaneous  aircraft  estimates  of 
erosion,  LWC,  and  cloud  thickness  quantities  were  calculated  by  means  of  a step- 
wise multiple  regression  program.  Since  this  program  assumes  linear  relation- 
ships, and  since  it  was  obvious  that  satellite  observed  brightnesses  and  tempera- 
tures were  not  linearly  related  to  erosion  parameters,  the  brightness  and  tempera- 
tures were  "linearized",  that  is,  transformed  into  numbers  which  are  linearly 
related  to  erosion  parameters.  It  was  also  noted  that  the  standard  deviations  and 
coefficients  of  variation  of  the  brightness  or  1R  temperature  were  not  significant 
predictors  of  erosion  parameters;  therefore,  only  the  averages  of  brightness,  1R 
temperature  and  IR  height  were  subsequently  linearized  and  used. 

Equations  of  the  following  forms  were  tested  to  develop  the  best  linearization: 

Y = A+B/X  , (6) 

Y = l/IA+BX)  , (7) 

Y = A (X)B  , (8) 

Y = A + B/(C-X)  . (9) 

In  general,  the  hyperbolic  form  of  Eq.  (9)  gave  the  best  fit  and  the  fraction,  1/(C'-X), 
was  used  to  linearize  X. 

As  an  illustration,  the  equation  found  to  linearize  brightness  13^  to  ESl^,  is; 

ESIp  = -2.554  + 506.253  /(288.  7 - T^)  . (10) 

Equation  (10)  and  the  data  points  from  which  it  was  derived  are  shown  in  Figure  17. 
The  correlation  coefficient  is  0.62  and  the  standard  error  of  estimate  is  1.61.  In 
this  case,  the  fraction,  1/(288.7-13^),  is  used  to  linearize  T3^  before  entry  into 
the  multiple  linear  regression  program.  The  program  then  derives  constants  for 
multiple  regression. 
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Figure  17.  Data  Points  and  Hyperbolic  Equation  Used  to 
Relate  Visible  Data  to  ESI 

P 

3.2  Dependent  Sample  Si7.e 

Equations  relating  satellite  parameters  to  erosion  parameters  were  first  de- 
veloped from  38  cases  with  simultaneous  aircraft  and  satellite  measurements  of 
clouds.  Due  to  pressure  for  early  results,  the  equations  were  applied  to  3 months 
of  satellite  data  over  Europe  and  Asia.  However,  four  additional  cases  having  very 
good  aircraft  data  became  available  later,  so  it  was  decided  to  include  these  as  part 
of  the  dependent  sample.  This  necessitated  recomputations  of  the  relation- 
ships and  using  the  new  relationships  to  replace  the  original  output  which  had 
been  delivered.  For  this  reason,  throughout  the  remainder  of  the  report, 
reference  will  be  made  to  equations  based  on  data  samples  of  different  size,  namely 
38,  41,  and  42  cases. 

The  most  desirable  prediction  equations  are  usually  found  bv  using  the  highest 
possible  number  of  cases  in  the  dependent  sample;  however,  in  the  present  sample 
the  last  case,  number  42,  was  extraordinarily  high  in  water  content  and  erosion 
parameters.  For  example,  when  ESI^  is  considered,  its  value  according  to  Table  4 
is  22.  20  as  compared  to  the  next  highest  value  of  7.44.  For  a sample  including 
this  extreme  case,  the  regression  program  will  tend  to  fit  an  equation  to  the  out- 
lying point,  at  the  expense  of  lower  values,  in  order  to  obtain  the  best  overall  fit. 
More  specifically,  when  regression  analysis  is  applied  first  to  41  cases  excluding 
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the  highest  case  and  then  to  all  42  cases  the  following  correlations  and  standard 
errors  of  estimate  were  found. 


Table  5.  Comparison  of  the  Correlation  Coefficients  (r) 
and  the  Standard  Errors  of  Estimates  for  ESI  Equations 


• * — — — 

4 1 Cases 

42  Cases 

r 

0.  7 85 

0.  915 

std  error 

1.  20 

1.  51 

At  first  glance  the  high  correlation  found  with  the  42  cases  is  very  attractive  and  the 

residual  error  for  the  42nd  case  is  near  zero.  On  the  other  hand,  the  standard 

error  of  estimate  for  the  42  case  equation  is  significantly  greater  than  that  for  t He 

41  case  equation  and  is  a serious  limitation  since  the  program  sponsors  specified 

that  maximum  accuracy  was  desired  in  the  2 to  8 range  of  ESI  . Similar  differences 

P 

between  the  41  and  42  case  equations  were  found  for  all  the  predictands  which  led 
to  the  adoption  of  the  41  case  sample  for  final  regression  analysis. 


2.2  Erosion  Liquid  Water  Content  and  ('.loud  Thickness  Ha-diet  ions 
From  brightness  and  IK  Data 

The  regression  equation  used  to  predict  ESI  from  brightness  and  If?  data, 
based  on  the  28  case  dependent  sample  is; 


ESI  = 


518 


14U 


IR  - 159.  3 27  l.  4 - B 


5.  84 


N 


(11) 


The  correlation  coefficient  is  0.79  and  the  standard  error  of  estimate  is  1.20. 

Similar  equations,  which  were  developed  to  predict  alt  the  erosion  parameters,  I VVC 

and  cloud  thickness  based  on  the  4 1 case  dependent  sample,  are  shown  in  Table  6. 

A plot  of  the  41  case  equation  for  certain  values  of  ESI  is  shown  in  Figure  18. 

P 

Aircraft  observed  ESIp  values  are  also  plotted  to  indicate  the  degree  of  scatter  and 
ranges  of  visible  anil  IR  covered  in  the  dependent  sample.  A dashed  line  shows  the 
approximate  range  of  satellite  data  from  a large  independent  sample  in  May  over 
Eurasia,  and  this  range  is  not  much  greater  than  the  range  of  aircraft  cases. 
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Table  6.  Equations,  Correlation  Coefficients,  and  Standard  Errors  of  Estimate 
for  Prediction  of  Erosion  Parameters,  LWC,  and  Cloud  Thickness  Based  on 
41  Cases 


Equation 

Correlation 

Coefficient 

Standard  Error 
of  Estimate 

ESI 

P 

746.249 

177.  500 

- 7.395 

0.785 

1.265 

(!ir  - 148.  800) 

(288. 7 00-ITpj) 

ESI 

P 

76. 981 

4.  932 

0.746 

1.359 

(18.  163  - Hm) 

ADI 

810.318 

81.354 

- 7. 317 

0.  690 

1.  634 

(IR-  - 150.400) 

(297.  600  - E"n) 

ADI 

71.371 

4.767 

0.  698 

1.  615 

(17.531  - H^l 

AR 

3138.756 

343. 004 

- 28.409 

0.  694 

6.  327 

(TIT  - 150.  800) 

(297.  900  - Bn) 

AR 

312.715 

19.985 

0.704 

6.  238 

(13.162  - 

AS 

284. 322 

83.7  10 

- 2. 110 

0.774 

0.  560 

OR  - 153.500) 

(288.900  - B^) 

AS 

34.  57  1 

1.  287 

0.741 

0.  594 

(18.  360  - Hy^) 

LWC 

476826. 657 

110445. 603 

- 2913.  73 

4 0.734 

362.  011 

(TR  - 7 1.460) 

(315.  680  - 

esit 

1065. 209 

+ 

266.  467 

- 9.  058 

0.793 

1.  287 

(TR'  - 129.  000) 

(297.400  -"Bn) 

AH 

209460. 983 

1200.  841 

- 156.  606 

0.  830 

1.  894 

(TRT  - 1097.  000) 

(355.  100  - ffN) 

Figure  18  like  all  the  erosion  equations,  shows  rapid  increase  in  ESI^  as 
brightness  approaches  its  maximum  and  as  IR  temperature  approaches  minimum 
values.  It  is  of  interest  to  note  the  behavior  of  the  ESI^  curves  for  273°,  the  melt- 
ing temperature  of  ice.  Precipitation  cases  with  cloud  top  temperatures  above 
273°K  will  have  few  or  no  ice  particles  so  that  high  ESI^  is  not  expected  in  mid- 
latitudes except  during  the  infrequent  cases  of  warm  rain.  The  curves  show  that 
zero  ESIp  is  found  with  a brightness  of  160  and  that  at  the  maximum  brightness 
found  within  the  data  sample,  namely  240.  the  ESI^  is  2.  Extrapolation  to  bright- 
nesses beyond  the  dependent  data  range  is  reasonable  since  very  bright  clouds  just 
above  freezing  may  represent  warm  rain  in  the  tropics  with  significant  ESl^  values. 
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Figure  18.  Plot  of  Equation  Used  to  Predict  ESIp  from  and 
TTTV  Based  on  41  Cases.  Dependent  data  are  sho\Vn 

3.4  Erosion.  Liquid  Water  Content,  and  Cloud  Thickness  Predictions 
From  IK  Data  Alone 

For  the  6-month  period,  November  through  April,  the  visible  brightness  over 
Eurasia  may  be  partly  a measure  of  the  reflection  from  snow  cover,  or,  in  mid- 
winter, it  may  be  unavailable  over  northern  portions  due  to  darkness.  These  limi- 
tations forced  us  to  consider  special  approaches  for  winter  conditions.  First,  a 
brief  experiment  was  run  to  determine  the  feasibility  of  using  brightness  for  at  least 
part  of  the  6-month  period.  Minimum  reflection  values  were  determined  for  6-day 
running  periods.  These  values  were  assumed  to  represent  ground  reflectance  and 
were  subtracted  from  the  daily  observations.  As  a result,  more  reasonable  erosion 
values  were  calculated  over  clear  areas  but  the  effect  of  decreasing  the  brightness 
of  heavy  clouds  yielded  ESI^  values  that  were  low.  It  was  also  found  that  over  the 
frozen  sea  or  tundra  snow  cover  alone  was  often  brighter  than  snow  cover  plus  thin 
cloud,  an  observation  not  understood  but  possibly  related  to  the  reflective  properties 
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of  snow  at  low  sun  angles.  In  view  of  these  complications,  it  was  decided  not  to  use 

the  brightness  in  any  form  over  Eurasia  during  the  6-month  period. 

Additional  calculations  revealed  that  the  use  of  IR  temperature  or  IR  height 

(11=-)  alone  gave  erosion  predictions  of  almost  the  same  accuracy  as  t,  - use  of  both 
IR 

IR  and  brightness  in  the  41  case  dependent  sample.  This  sample  was  from  a rela- 
tively homogeneous  area  over  winter  and  spring  where  temperature -height  profiles 
did  not  vary  appreciably.  However,  in  view  of  the  much  greater  variation  of 
temperature -height  profiles  over  Eurasia  during  the  6-month  period,  the  use  of 
H-j^-  appeared  to  be  more  desirable  than  IR  temperature  alone.  Therefore, 
plus  TR"  predictors  were  used  for  the  months  of  May  to  October  1973,  while  the 
predictor  was  used  for  the  months  of  April,  November,  December  1973,  and 
January  to  March  1974.  Equations  using  alone  to  predict  all  the  erosion  param- 
eters, LWC,  and  cloud  thickness  from  the  41  case  dependent  sample  are  given  in 
Table  6.  ITT  temperatures  were  converted  to  Hj^-by  means  of  average  monthly 

temperature -height  profiles.  These  profiles  were  developed  by  Spiegler  and 
33 

Fowler  from  a 4 to  5 year  data  bank.  A profile  was  made  for  every  32  original 
satellite  grid  points  (the  NMC  grid).  In  practice  each  5X5  data  box  used  the  nearest 
profile  which  was  never  more  than  22.  6 original  grid  points  distant.  At  latitude 
65°  this  is  about  27  1 km.  Larger  errors  in  the  converted  value  are  likely  to 

occur  due  to  a day's  departure  from  the  average.  Examination  of  standard  devit 
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tions  at  standard  pressure  levels  as  given  by  Crutcher  and  Meserve  by  month  and 
area  show  a maximum  of  about  4°  in  April  at  700  mb  (3  km)  over  central  Eurasia. 

In  the  Moscow  area  the  maximum  is  about  3°  at  700  mb  (3  km)  and  at  500  mb  (5.  5 km), 
also  in  the  spring  months.  These  temperature  variations  are  equivalent  to  about  0.7 
and  0.  5 km  respectively,  and  when  placed  in  the  erosion  equations  give  the  range 
of  uncertainty  listed  in  Table  7. 


Table  7.  Erosion  Parameter  Uncertainties  Corresponding  to  Maximum 
Monthly  Standard  Deviations  of  Temperature  at  Two  Altitudes 


at  3.  0 km 

at  5.  5 km 

± 0.  1 

± 0.  2 

± 0.  2 

±0.3 

± 1.0 

±1.0 

± 0.  1 

± 0.  1 

33.  Spiegler,  D.  B. , and  Fowler,  M.  G.  (1972)  Four  Dimensional  World  Wide 

Atmospheric  Models  (surface -25  km  altitude)  NASA  Crt-iS0ft2,  Environmental 
Res.  and  Tech. , Inc.,  Concord,  Massachusetts  01742,  65  pp. 

34.  Crutcher,  H.  L, , and  Meserve,  J,  M.  (1970)  Selected  Level  Heights,  Tempera  - 

ture  and  Dew  Points  for  the  Northern  Hemisphere.  NAV AlK  hO-lC-52. 
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It  is  therefore  concluded  that  use  of  average  monthly  temperature-height  profiles 
introduces  relatively  small  errors  while  saving  significant  time  and  expense. 

Determination  of  heights  above  the  ground  from  the  0 to  15  km  profiles  were 
made  according  to  the  following  rules: 

(1)  First,  determine  tropopause  height. 

(a)  If  change  of  slope  AH  AT  becomes  - 0.25  km  °K  * from 

1 km  level  to  the  next  going  upward,  this  level  is  assigned 
the  tropopause  height. 

(b)  If  not.  then  level  of  lowest  temperature  is  designated  as 
the  tropopause. 

(21  Height  of  IR  temperature  is  then  determined  progressing  downward  from 
the  tropopause  level. 

(a)  If  IR  temperature  exceeds  any  in  profile,  assign  to  level  of  highest 
temperature  in  profile. 

(b)  If  IR  temperature  is  lower  than  any  in  profile,  assign  to  level  of 
of  lowest  temperature. 

Further  tests  revealed  a bias  when  erosion  parameters  were  calculated  over  a 
wide  range  of  conditions  from  both  the  plus  TTT  and  from  the  H|^- alone  equations 

given  in  Table  6.  The  H— alone  equations  gave  higher  ESI  in  heavy  weather. 

IK  P 

Therefore,  in  order  to  maxe  the  estimates  over  winter  and  summer  as  homogeneous 
as  possible,  erosion  parameters  were  calculated  from  both  and  TTT and  from 
1 1^- alone  for  a large  sample  and  compared.  The  sample  was  chosen  to  contain  a 
minimum  of  zero  erosion  values.  A comparison  of  the  two  equations  was  made  for 
each  erosion  parameter  over  a critical  range  of  values  for  that  parameter.  In  com- 
paring ESIp  from  the  two  equations  for  example,  cases  with  ESl^  between  0 and  5 
according  to  the  T5^,  plus  TTT  equation  were  used.  This  reduced  the  sample  size 
somewhat.  The  comparison  consisted  of  developing  a new  linear  regression  equa- 
tion between  the  two  sets  of  ESI^  values  from  which  an  intercept  and  slope  could  be 
obtained. 

The  equations  for  converting  erosion  parameters  derived  from  Hy^-onlv  to  TR 
and  equations  are  listed  in  Table  8.  Cases  outside  the  limits  given  were  ex- 
cluded. Correlations,  standard  errors  of  estimate,  and  numbers  of  cases  in  the 
sample  are  also  listed.  In  all  cases,  the  erosion  parameters  from  the  11—  variable 
alone  must  be  reduced  for  best  agreement  with  the  and  TTT  equations.  The  most 
probable  cause  is  overemphasis  to  cirrus  cases  in  the  11^  only  estimates.  TTT  and 
histograms  over  large  areas  show  numerous  cold  but  medium  bright  cases  of 
cirrus  which  probably  are  underrepresented  in  the  4 1 case  sample  used  to  form  the 
original  equations.  The  equations  in  Table  6 were  combined  with  those  in  Table  8 
to  provide  the  final  equations  which  are  shown  in  Table  A plot  of  final  ESI 
equation  for  H^-  alone  is  shown  in  Figure  18.  The  figure  shows  that  ESI^  values 
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of  zero  occur  for  systems  having  values  up  to  2.  6 km.  In  the  area  of  lowest 
mid-winter  Furasian  tropopause  heights,  where  maximum  is  about  8 km,  the 
ESlp  will  not  exceed  1.8.  Southern  tropopause  heights  of  15  km  will  give  an  ESI^ 
of  about  12.5  assuming  that  dense  clouds  are  measured  at  that  level. 


Table  8.  Conversions  from  Equations  using  H*»j-Only  to  Equivalent  Equations  of  ITT 
and  Bn  IH 


Equation 

Upper 
Lim  it 

Correlation 

Coefficient 

Standard 
Error  of 
Estimate 

% 

Number 
of  Cases 

ESi  = 0.  679(ESI  ) + 0.  045 

P P 

mVRN  Hjg- 

S\ 

8.0 

0.  892 

0.  582 

26,  850 

ADI  = 0.  600( ADI  ) + 0.  105 

"nr 

/\ 

11.0 

0.  895 

0.  584 

26, 885 

AR  = 0.  640(AR  ) + 0.  300 

Hnr 

/N 

50.  0 

0.  907 

2.  189 

26,  885 

AS  = 0.  920(AS  ) - 0.274 

nr+BN  FInr 

4.  0 

0.  893 

0.  355 

26,  816 

For  7 and  8 May  1973  Total  data  points  = 27,  626 

Satellite  data  missing  = 741 

Satellite  data  O.K.  = 26,885 


Table  9.  Final  Equations  to  Predict  Erosion 

Parameters  from  H— Alone,  Based  on  41  Cases 
IR 

ESIp  = (lB.'lS?  - H^)  ' 3,304 


42.  823 

(T775FT-  Hjy) 

200. 138 

(18.162  - Hj^-) 

31.  805 

(TOTBTT  - H^-) 


- 2.755 

- 12.490 

- 1.458 
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Figure  19.  Plot  of  Final  Equation  Used 


to  Predict  ESI  from 
4 1 cases 


H 


ITT 


Based  on 


3.5  Probability  Prediction* 

A probability  model  by  Gringorten5  was  applied  to  the  data  in  Figure  18  in  order 
to  establish  the  probabilities  that  a threshold  of  ESIp  would  be  exceeded  as  a func- 
tion of  cloud  temperature.  TfT.  and  brightness,  BN.  Figure  20  has  the  results  for 

an  ESI  threshold  of  5.  Probabilities,  based  on  the  38  case  sample,  are  as  high  as 
p 

90  percent  for  the  coldest  and  brightest  clouds  observed. 


> 


PROBABILITY  ESI  PRECIP  i 5 


IR  *K 


Figure  20.  Probabilities  of  Exceeding 
an  ESI  of  5 Related  to  Satellite 
Measurements 
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4.  APPLICATION— CLIMATOLOGICAL  AN1)  REAL  TIME 
EROSION  PARAMETERS 


4.1  Program  Schedule 

In  response  to  a letter  from  SAMSO,  dated  7 May  197  6,  in  which  various  cli- 
matologies of  erosion  parameters  were  requested,  the  program  schedule  in  Table  10 
was  drawn  up  and  approved  by  23  June  197  6. 


Table  10.  Program  Schedule  for  Satellite  Estimates  of  Erosion  Parameters 


T ask 

Date  Due 

Date  Mailed 

i. 

Method  Verification 
(as  developed) 

1 July  197  6 to 

1 October  1977 

2. 

ESI  maps 

p r 

May  to  October  1973  (6  mos) 

April  1973,  November  to 
December  1973 
(remaining  6 mos) 

January  to  March  1974 

15  December  1976 

1 October  1977 

10  December  197  6 

Cancelled 

3. 

Target  statistics 

May,  July,  October  1973  (3  mos) 
Remainder  (9  mos) 

1 October  197  6 

1 October  1977 

18  October  1976 

13  December  19*6 

30  October  1977 

4. 

ESI  - test  areas 

■ 

Wallops  & Kwajalein 

May  to  October  1973  (6  mos) 
Remainder  (6  mos) 

1 October  197  6 

1 March  1977 

10  October  197  6 

8 March  1977 

5. 

ADI,  AR.  AS 

correlations 

Test 

Remainder  (12  mos) 

1 March  1977 

1 October  1977 

Verbal  ok  21  April  1977 

3 October  1977 

3 mos  only 


Over  the  work  period  changes  occurred,  new  products  were  generated  and  some 
of  the  original  products  requested  were  simplified  but  in  general  the  program  re- 
mained intact.  Completion  dates  for  the  various  parts  of  the  program  are  also  shown 
in  Table  10. 

4.2  Grid  Siw 

Most  of  the  dependent  cases  from  which  the  regression  equations  were  developed 
used  satellite  data  averaged  over  50  to  70  samples  usually  in  parallelogram  shaped 
areas.  These  are  roughly  equal  in  size  to  square  areas  of  8 X 8 data  points.  It  was 
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decided  that  useful  applications  of  the  erosion  equations  could  be  obtained  from 
5X5  square  areas,  although  the  significance  of  equations  when  applied  to  smaller 
areas  was  unknown  at  the  time.  Therefore,  unless  otherwise  noted,  all  erosion 
parameters  derived  from  the  NOAA  satellite  archive  are  for  5X5  squares  which 
range  in  actual  size  from  33  X 33  km  at  the  equator  to  61  X 61  km  at  70°N.  Further- 
more, it  was  arbitrarily  decided  that  a valid  5X5  computation  required  data  from 
at  least  15  of  the  25  possible  data  points. 

4.3  Km ironmenlal  Severity  Index-Precipitation  (K.SI^)  Maps 

ESI  maps  and  target  statistics  are  described  in  a separate  annex  to  this  re- 
port. The  annex  is  available  to  qualified  agencies  upon  request  to  AFGL/LY, 
Hanscom  AFR,  MA  0173  1. 

4.4  K.SIp  for  tlie  Knajaleiii  ami  Wallops  Test  Kanges 

ESI  values  over  5X5  grid  areas  were  computed  from  the  38  case  plus  ITT 
equation.  Brightness  was  used  year  around,  even  at  Wallops,  since  the  prime  areas 
of  interest  lay  over  the  sea.  Arrays  of  values  were  printed  out,  one  day  per  page, 
and  grid  overlays  were  supplied  for  positioning  the  data. 

A sample  printout  for  the  Wallops  area  with  superimposed  overlay  is  shown  in 
Figure  21.  The  geographical  grid  is  skewed  in  relation  to  the  satellite  grid  as 
discussed  in  Section  2.  1.  1.  Missing  data  are  indicated  by  999.  A similar  printout. 
Figure  22,  shows  ESI^  for  the  Kwajalein  area. 

In  the  Wallops  area,  satellite  data  were  frequently  missing  because  a satellite 
readout  station  is  located  at  that  site  and  interrupts  local  data  when  reading  out 
stored  orbits.  A count  of  useful  ESI^  arrays  at  both  test  ranges  is  given  in  Table  11. 

4.5  Kraltinit-  Application  of  ESIp  Erosion  Parameter 

Once  it  became  possible  at  AFGL  to  accurately  navigate  data  from  the  geo- 
synchronous satellites,  a program  was  written  to  normalize  the  brightness  and 
compute  ESlp  in  realtime.  All  grid  data  were  for  the  nominal  8 km  subpoint  reso- 
lution discussed  in  Section  2.  1.  2. 

An  8 X 8 array  of  ESI^  values  could  be  presented  on  the  McIDAS  CRT  screen 
3 min  after  the  satellite  scanned  the  clouds.  A 50  X 50  array  could  be  printed  out 
and  positioned  under  a grid  in  12  to  13  minutes. 


' 
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WALLOPS  IS 


WALLOPS  OVERLAY 


Figure  21.  ESIp  in 
NOAA-2  Satellite 


the  Vicinity  of  Wallops  Island  1 April  1973.  fror  the 


The  ESIp  data  and  cloud  motions  were  frequently  relayed  by  phone  to  the  V i 
Weather  Team  at  Wallops  in  the  spring  of  1977.  Generally,  positions  were  . 
in  terms  of  azimuth  and  distance  from  Wallops.  The  ESI  estimates  were  use 
approximate  erosion  in  areas  beyond  reach  of  the  SPANDAR  radar.  In  additi  - 
for  the  first  time,  comparisons  of  ESI^  values  obtained  by  radar  and  sate  'it* 
come  available  in  areas  closer  to  the  SPANDAR.  Unfortunately,  onlv  a few 
titative  comparisons  were  obtained  but  the  team  did  become  accustomed  to  th. 
advantages  and  limitations  of  the  satellite  estimates. 
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KWAJAlflN 


?I!KW  V APR  rs 


KWAJALEIN  OVERLAY 


Ktgure  22.  KSl„  in  the  Vicinity  of  Kwajalein  Atoll  9 April  1972  from 
the  NO  A A -2  SatMHte 

( leosynchronous  satellite  estimates  of  KSI  over  and  upwind  of  the  Kwajalein 
area  were  first  thought  to  he  useless  due  to  a -zenith  angle  of  <? 5 ° to  the  nearest 
satellite  at  0*N,  1 9 !)'  W . However,  KSI  values  were  calculated  from  the  eastern 
satellite  for  a small  area  which  duplicated  the  viewing  and  sun  angles  lo  Kwajalein 
from  the  western  satellite  at  the  time  of  maximum  interest.  These  were  com- 
pared with  values  calculated  15  min  earlier  for  the  same  area  from  the  western 

satellite  which  looked  down  at  only  12*  from  the  vertical.  The  KSI  values  were 

P 

surprisingly  close,  with  maximum  values  in  individual  clouds  differing  by  only  a 
few  percent.  Naturally.  1 lie  scales  were  different  and  narrow  K-W  variations  in 
the  clouds  were  not  evident  when  viewed  at  the  high  zenith  angle.  Kxamples  are 
presented  In  Section  !>,  5 of  this  report.  After  this  comparison  a procedure  for 
relaying  KSI  patterns  and  cloud  motions  for  realtime  use  at  Kwajalein  began. 


Table  11.  Frequency  of  Satellite  Coverage  at  Test  Hanges 
During  One  Year 


Wallops  Days 

With  Data 

Kwajalein  Days 
With  Data' 

April 

1973 

18 

28 

May 

1073 

10 

31 

.lune 

1073 

25 

20 

duly 

107  3 

24 

20 

August 

1073 

22 

3 1 

September 

107  3 

24 

20 

October 

1073 

26 

30 

Novemhe  r 

1073 

24 

25 

December 

107  3 

22 

30 

January 

107  4 

23 

30 

February 

1074 

25 

27 

March 

107  4 

16 

26 

Total  268 

345 

% Possible  7 3 

05 

Since  the  western  geosynchronous  satellite  does  not  oscillate  appreciably  from  Die 
eipiator  in  its  orbit  it  is  possible  to  use  master  grids  in  equatorial  areas.  MelUAS 
programs  calculate  latitudes  and  longitudes  of  the  center  and  corner  (mints  of  arrays 
of  satellite  data.  These  indices  are  then  used  for  positioning  the  overlay  grids.  The 
most  useful  data  were  relayed  near  20110/  or  about  11  hr  before  the  usual  launch 
window  time.  Data  for  upwind  areas  beyond  radar  range  proved  useful.  Once 
again  comparisons  between  radar  and  satellite  derived  FSl  were  obtained  within 
radar  ranges.  A program  of  taping  the  data  during  the  eastern  night  was  developed 
to  facilitate  subsequent  analysis  and  comparisons.  A sample  printout  of  l-;SI  values 
and  grid  overlay  is  shown  in  Figure  2.1. 
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statistically  l\v  dividing  the  sample  into  two  parts.  This  approach  allowed  equations 
to  lie  checked  against  independent  data  and  predicts  how  well  the  equations  w ill  work 
when  applied  to  other  sets  of  independent  data  such  as  the  satellite  data  over  target 
areas.  Second,  satellite  estimates  were  compared  with  a limited  sample  of  ESI^ 
estimates  from  radar  at  temperate  and  tropical  latitudes.  This  approach  is  unique 
in  that  the  radar  ESI  estimates  are  instantaneous  and  synchronous  with  the  satel- 
lite. However,  individual  samples  were  for  horizontal  areas  whose  size  was  con- 
siderably less  than  the  areas  viewed  by  the  satellite  sensors.  Third,  satellite 
estimates  were  compared  with  a limited  sample  of  ESI^  estimates  from  aircraft 
soundings  from  eonvei  Live  systems  near  Kwajalein.  In  this  ease,  the  technique  for 
obtaining  the  mass  profile  from  aircraft  was  the  same  as  in  the  41-ease  sample, 
but  the  comparisons  yield  useful  information  on  the  value  of  ESI  estimates  in  highly 
convective  areas  which  were  not  included  in  the  41-case  sample.  Fourth,  another 
comparison  shows  the  effects  of  slant  viewing  from  a satellite  on  ESI  values  near 
Kwajalein.  Fifth,  a large  sample  of  satellite  estimates  were  compared  with  near 
synchronous  ESI^  estimates  from  the  A F(1 1 , 11  model  at  11  stations  in  Europe  and 
Asia.  This  comparison  gives  no  information  on  the  absolute  accuracy  of  either 
method.  It  merely  compares  two  different  techniques  for  estimating  ESI  . 

There  are  limitations  to  verifications  which  are  important  to  remember.  The 
most  important  limitation  is  that  no  accepted  standard  exists  for  measuring  hvdro- 
meteor  distributions  over  substantial  volumes,  l'lic  satellite  correlation  studv  has 
used  procedures  developed  very  recently  by  AFtll.  to  support  reentry  system  ero- 
sion tests  and  this  technology  is  still  undergoing  refinement.  Consequently,  the 
extent  of  uncertainty  in  the  hydrometeor  measurements  described  is  not  well  known. 
Another  limitation  is  the  fact  that  the  satellite  measurements  of  clouds  have  not 
been  compared  directly  to  measurements  of  nosecone  erosion  in  clouds  but  rather 
to  parameterized  hvdrometeor  profiles  (ESI  , Al)l)  or  calculated  erosion  effects 

p 

(AS,  All).  How  the  uncertainties  in  modelling  erosion  interact  with  uncertainties 
in  satellite  and  aircraft  measurements  is  not  known. 


5.2  Statistical  Approach 

A common  approach  in  verification  of  empirical  relationships  is  to  divide  avail- 
able data  into  two  mutually  exclusive  sets  named  the  dependent  set  and  the  indepen- 
dent set.  A relationship  is  developed  from  the  dependent  set  and  then  applied  to  the 
independent  set  as  a test.  The  errors  generated  in  the  test  predict  the  expected 
accuracy  in  future  use  of  the  relationship  with  new  sets  of  data. 

A slight  departure  from  this  approach  was  necessary  due  to  the  limited  number 
of  dedicated  and  expensive  aircraft  flights  available  to  our  study.  Only  4 1 eases  were 
available  when  final  applications  had  to  be  made.  It  was  decided  that  the  highest 
payoff  would  result  by  using  all  the  cases  to  derive  relationships  rather  than  divide 


57 


1 


the  data  and  risk  insignificant  relationships.  Applications  proceeded  with  rela- 
tionships based  on  41  cases  as  discussed  in  Section  3. 

A later  experiment,  however,  revealed  that  similar  relationships  could  be 
derived  from  fewer  cases  so  t hat  the41  case  sample  could  be  divided  and  relation- 
ships could  be  checked.  One  of  the  41  cases  having  a very  low  ESl^  was  selected 
at  random  and  removed  to  give  an  even  number  of  cases.  The  cases  were  then 
ranked  in  terms  of  aircraft  10SI^  and  subdivided  into  4 groups  of  10  each  so  that 
prediction  c»rrors  from  high  and  low  ESI^  cases  could  be  studied  separately.  The 

Hoot  Mean  Square  (rmsl  residuals  in  Table  12  in  the  column  marked  "dependent" 

❖ 

are  the  results  when  the  38-case  ESI^  equation  is  used  on  the  40-case  sample. 
The  residuals  are  simply  the  differences  between  the  aircraft  observations  of  ESI 
and  the  predictions  of  ESI^  from  the  equation.  The  residuals  increase  as  ESI^ 
increases;  however,  the  quotient  of  the  residual  divided  by  the  ESl^  decreases  as 
KSIp  increases.  In  the  important  range  of  ESl^  from  2 to  8,  residuals  are  less 
than  the  observed  ESl^.  The  residuals  for  API,  AS,  and  AH  exhibited  similar 
tendencies. 


Table  12.  Hoot- Mean-Square  (rms)  Residuals  of  ESlp  in 
Dependent  (40  case  sample)  and  Independent  (2  X 20  case 
samples) 


Ranked 

Samples 

Range  of 
Aircraft  ESlp 

rms  Residuals 

Dependent 

Independent 

Nos.  1 to  10 

0.0  to  0.06 

0.  64 

0.  85 

Nos.  1 1 to  20 

0.  07  to  0.  63 

0.  66 

0.  76 

Nos.  21  to  30 

0.  74  to  2.  64 

1. 04 

1.  24 

Nos.  31  to  40 

2.  68  to  7.  44 

2.  24 

3.  66 

Nos.  1 to  40 

0.  0 to  7.  44 

1.32 

2.  01 

The  residuals  in  Table  12  column  marked  "independent"  are  results  when  the 
40-case  sample  is  randomly  divided  in  half,  a prediction  equation  derived  for  each 
half  separately  and  then  used  to  predict  ESI  values  for  cases  in  the  other  half- 
sample. The  independent  residuals  are  higher  than  the  dependent  residuals,  as 
expected;  however,  since  they  arc  not  appreciably  higher  and  since  all  the  predic- 
tion equations  are  similar  we  expect  that  the  extensive  applications  reviewed  in 

The  38-case  equation  is  reasonably  close  to  a best  fit  to  the  40-case  sample  since 
the  residuals  of  the  last  two  cases  from  the  predictions  of  the  38-case  equation  are 
not  large.  Consequently,  no  special  equation  was  derived  for  the  40-case  sample. 
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Section  4 have  errors  which  do  not  exceed  the  independent  residuals  given  in 
Table  12.  The  results  of  Table  12  further  imply  that  meaningful  relationships  could 
have  been  developed  from  a sample  of  cases  half  the  size  of  the  one  available  to  this 
study. 

5.3  Radar  vs  Satellite  KSlp  Values 

5.3.  1 SPANDAR  RADAR  vs  GOES-EAST  SATELLITE- 
24  FEBRUARY  1977 

RHI  data  from  the  SPANDAR  radar  at  Wallops  Island  were  collected  from  19 

to  23Z.  Scans  not  more  than  3 min  from  satellite  time  were  then  processed  to 

give  ESlp  values  at  varying  distances  between  24  and  96  km  from  the  radar.  Levels 

of  particle  types  and  their  transition  zones  were  assigned  by  Dr.  R.  M.  Cunningham 

at  the  time.  These  were  based  on  air  temperature  and  aircraft  observation  of  the 

particles.  Equations  to  convert  radar  reflectivities  to  mass  densities  for  various 

24 

particle  types  have  been  summarized  by  Plank.  During  the  first  part  of  the 
period,  echoes  were  found  between  about  3 and  9 km  altitude  with  radar  ESI^  of 
0.  1 to  2.  2.  Near  the  middle  of  the  period,  precipitation  reached  the  ground  and  a 
maximum  ESlp  of  12.4  by  radar  was  computed.  Shortly  before  the  trend  to  higher 
ESlp  daylight  waned  and  the  comparisons  had  to  be  restricted  to  satellite  ESl^  by 
alone. 

The  comparison  is  summarized  in  Table  13.  The  column  marked  "n"  refers 

to  number  of  samples.  During  the  first  portion  of  the  period,  1904  to  2106Z,  tops 

by  satellite  averaged  0.  5 km  above  the  radar  and  ESI's  both  by  4 HF  and  lly^ 

alone  averaged  about  0.7  above  those  bv  radar.  The  RMSE  was  1.3  to  1.4  ESI^. 

During  the  second  period,  tops  by  satellite  averaged  1.  1 km  lower  than  the  first 

period.  Aircraft  reports,  although  not  exactly  over  the  sample  area,  also  showed 

a decline  in  top  altitude;  however,  the  average  radar  tops  rose  0.  1 km  from  the 

first  to  the  second  period.  This  may  be  attributed  to  echoes  received  by  the  lower 

side  lobe  of  the  antenna  pattern  being  misidentified  as  higher  echoes  in  the  main 

lobe.  KSlp  values  by  radar  were  much  higher  than  those  by  satellite  during  the 

second  period.  Examination  of  the  mass  profiles  shows  that  this  is  only  partially 

due  to  the  side  lobe  effect  because  very  high  water  content  was  found  below  3 km, 

the  freezing  level.  The  cause  for  these  large  differences  is  not  known  nor  do  we 

know  how  often  this  occurs  since  the  sample  is  small. 

Availability  of  visible  data  through  the  second  period  might  very  well  have 

raised  the  satellite  estimates  of  ESI  but  they  would  still  be  less  than  ESI  esti- 
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mates  bv  radar. 
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5.  3.  2 TRADEX  RADAR  vs  GOES-WEST  SATELLITE  — 
8 JULY  1977 


This  comparison  was  made  near  Kwajalein  Atoll  through  clouds  at  the  northern 
edge  of  an  Inter  Tropical  Convergence  Zone  cloud  system.  ESI^  values  were  com- 
puted along  an  RHI  obtained  at  0210Z  on  the  120°  radial  from  the  TRADEX  radar. 

Mr.  C.  N.  Touart  read  reflectivities  by  eye  from  a photograph  of  the  color-con- 
toured RHI  display.  Profiles  of  mass  at  1 km  intervals  were  calculated  and  ESI^ 
values  obtained.  Figure  24  compares  these  values  with  ESI^  derived  from 
GOES-W,  which  scanned  the  area  at  0223Z.  The  41  case  equation  with  variables 
"Bpj  and  TR  was  used  to  calculate  ESI^.  The  GOES  navigation  at  this  time  was  good 
to  about  54  km.  The  variation  of  radar  values  along  the  radial  and  the  variation  of 
the  satellite  contours  match  well.  The  maximum  of  4.  1 on  the  RHI  was  in  a shower 
area  while  the  minimum  of  1.  2 occurred  with  cloud  between  8 and  12  km.  The  maxi- 
mum of  3.  8,  found  furthest  from  the  radar,  had  a radar  horizon  near  8 km,  but 
showers  probably  fell  to  the  ocean.  Adjustments  were  not  applied  for  the  missing 
layer  from  the  radar  horizon  to  the  sea,  but  if  they  had  a slight  increase  would 
occur  at  the  maximum  of  4.  1 and  probably  a significant  increase  at  the  maximum 
of  3.  8. 


It  should  be  remembered,  as  noted  in  Section  2,  that  under  these  satellite  view- 
ing conditions  a cloud  12  km  high  is  positioned  about  24  km  too  far  to  the  west.  This 
case  improves  the  comparison  by  placing  the  higher  radar  values  even  closer  to  the 
higher  satellite  values.  The  comparison  is  considered  good,  especially  in  view  of 
the  uncertainties  noted. 


5.4  Aircraft  v«  Satellite  F.SIp  Values  Near  Kwajalein 

5.4.  1 THE  8 AUGUST  1976  CASE 

This  was  a case  of  a small  isolated  shower,  with  anvil,  but  without  secondary 
clouds.  The  observational  data  were  excellent.  The  Kwajalein  WSR-57  radar  data, 
which  showed  other  isolated  showers,  could  be  used  to  precisely  locate  the  satel- 
lite data.  In  addition  to  good  1-D  PMS  aircraft  data  during  ascent  and  descent, 
excellent  time  lapse  and  still  photography  made  it  possible  to  develop  detailed  plots 
of  cloud  positions  and  shapes.  The  cloud  underwent  a minimum  of  change  during 
the  soundings.  It  drifted  slowly,  shear  was  weak,  and  mass  measurements  were 
essentially  the  same  during  ascent  and  descent. 


I 
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Figure  24.  ESIp  (in  tenths)  Derived  from  ClOES-W.  Each  digit  represents  a 
rectangular  area  8 km  N-S  by  15  km  W-E  (left  - right).  ESIp  (in  tenths)  cal- 
culated from  the  TRADEX  radar  is  shown  along  the  radar  radial 

The  cloud  appeared  to  be  "pumping",  that  is,  one  cell  after  another  grew  ver- 
tically and  spread  at  upper  levels.  The  spreading  levels  are  shown  in  Figure  25, 
a photograph  from  the  8.  0 km  level.  The  full  vertical  extent  of  the  cloud  is  de- 
picted in  cross  section  along  the  aircraft  track  in  Figure  26.  Precipitation  nearly 
evaporated  by  the  time  it  reached  2 or  3 km  altitude  only  to  be  replaced  by  another 
growing  cell.  Principal  aircraft  penetrations  were  in  the  4.5  to  6.0  km,  9.2,  and 
12.  8 km  levels.  Ry  interpolation  between  levels  and  the  use  of  visual  clues,  mass 
profiles  were  generated  for  the  core  of  the  cloud  and  several  distances  out  from  the 
core.  These  profiles  are  also  shown  in  Figure  26.  From  these  ESI^  was  calculated. 
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Figure  25.  Cumulonimbus  Anvil  Showing  Two  Levels  of 
Spreading  Below  the  Top.  Taken  fiom  8.  0 km  MSL, 
about  7 km  distant.  1 -D  probes  show  on  the  front  of  the 
aircraft 


Figure  26.  Left:  Mass 
Profiles,  at  the  Storm 
Core,  3 and  6 km  NE  and 
SW  from  the  Core  Along 
the  Aircraft  Track, 

Right:  Cross  Section  of 
the  Storm  Along  the 
Aircraft  Track 
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The  satellite  data  mesh,  radar  echoes,  aircraft  track  and  ESIp  values  were 
then  mapped  and  superimposed  as  shown  in  Figure  27.  Satellite  and  TR”  data 
from  the  NOAA  4 satellite  for  each  square  were  used  to  compute  ESIp  as  shown 
using  the  41  case  equation. 


Figure  27.  WSR -57  Radar 
Echo  (stippled)  of  the 
Storm  at  Satellite  Pass 
Time,  ESI„  Values  Along 
the  Aircraft  Track  and 
Estimated  ESI  Contours. 
Satellite  ESIp  in  italics 
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In  the  SW  the  anvil  lav  over  clear  air;  to  the  NE  small  cumulus  were  found 

under  the  anvil.  Agreement  between  the  satellite  and  aircraft  derived  ESIp  values 

is  good  over  both  areas  of  the  anvil.  The  core  which  had  an  ESI  of  15.  2 was  esti- 

P 

mated  to  be  only  2 km  wide  as  opposed  to  the  7. 5 X 7.  5 km  satellite  archive  mesh. 
For  this  viewing  geometry,  the  IR  sensor  resolution  is  about  8.2  km  and  the  visible 
resolution  is  half  that  size  so  these  could  not  be  expected  to  characterize  the  core 
of  the  storm.  The  satellite  archive  square  which  contained  the  storm  core  had  an 


KSIp  of  2.  4 while  the  average  of  the  aircraft  ESI^  contours  over  that  square  was  about 

5.  5 ESI  . Thus  the  average  ESI  was  2 times  the  satellite  estimate  while  the  core 
P P 

ESI  was  6 times  the  satellite  estimate.  Table  14  summarizes  the  comparisons. 

P 

Satellite  estimates  for  adjacent  squares  were  weighted  for  comparison  to  aircraft 
estimates  in  between  squares. 


'Fable  14.  ESIp  Determined  bv  Aircraft  and  bv  Satellite  in 
Relation  to  the  Core  of  a Small  Isolated  Cumulonimbus  Cloud 


Location 

Satellite 

Shower  core 

2.  4 

3 km  NE,  edge  of  precip. 

2.  0 

6 km  NE,  some  ci,  cu  below 

0 

3 km  SW  in  precip. 

1.7 

6 km  SW  anvil  only 

1.  1 

12  km  SW  anvil  only 

0.  5 

Brightest  and  coldest 
(combined  adjacent  sat. 
squares) 

2.9 

Aircraft 


5.4.2  THE  12  AUGUST  1976  CASE 

In  this  case  the  cloud  system  was  more  complicated  than  the  8 August  197  6 
cloud  system.  The  cloud  sampled  was  a dying  Cb  whose  base  was  about  6 km  alti- 
tude. It  was  located  on  the  northern  edge  of  a cluster  of  Cb's  of  which  two  were 
active  with  heavy  showers.  The  Lear  aircraft  PMS-lDand  photography  were  good. 
A happenstance  photo  of  the  aircraft  radar  display  proved  vital  when  compared  with 
the  WSR-57  radar  data  to  verify  a cloud  position. 

The  mass  profile  shown  in  Figure  28  was  a smoothed  fit  of  observed  densities 
to  give  what  was  considered  to  be  the  most  representative  profile  through  the  sound- 
ing area.  From  this  an  ESI  of  4.7  was  calculated. 

P 

A printout  of  ESI^  values,  their  contours,  and  WSR-57  radar  echoes  is  shown 
in  Figure  29.  The  ESI^  values  were  calculated  from  NOAA  4 satellite  data  using 
the  41  case  equation.  The  racetrack  sounding  area  of  the  aircraft  is  enclosed  in 
a box-like  area  near  ION,  168E.  Average  ESI^  within  this  box  is  4.  6,  a value 
very  close  to  that  calculated  from  the  aircraft  data. 
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Figure  28.  Mass  Profile  and  Observer  Notes  Through 
Dying  Cumulonimbus  (Cb)  Near  Kwa.jalein  on  12  August 
197  6.  Points  are  PMS-1D  system  measurements 


Comparison  of  the  WSR-57  radar  echoes  based  on  beam  scans  at  0“  elevation 

and  the  ESI  contours  is  of  interest.  In  the  sample  area  no  echo  was  seen.  The  lack 
P 

of  echo  is  reasonable  since  the  aircraft  observed  no  low  precipitation  with  the  com- 
plete loss  of  virga  at  4.  8 km.  The  echo  near  9.  5N  and  167.  5E  conforms  quite  well 

with  heavy  satellite  FSI  except  to  the  north  where  it  covered  some  ESI  values  of 
P P 

zero.  In  that  area  was  as  low  as  38  and  the  1R  temperature  as  high  as  286I,K. 

At  this  range  of  110  km  from  the  radar,  the  beam  would  not  see  the  lowest  altitudes 

and  therefore  could  not  detect  rain  from  low  ctimulus.  More  likely,  the  echo  was 

from  virga  with  little  or  no  cloud  below  3 km.  It  was  obviously  quite  transparent 

to  the  IR  and  not  very  reflective  in  the  visible.  Near  9.  ON,  168.  0E  a portion  of 
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another  echo  is  found  in  an  ESI  area  of  zero.  These  clouds  had  aH.  ' 

^ P N 

while  the  IK  was  as  low  at  27  1,,K.  The  echo  here,  only  40  km  from  the  i , 1 , 


most  likely  due  to  rain  from  cumulus  that  were  not  quite  bright  or  . old  • ■ i. 

give  an  ESI^  above  zero  according  to  the  general  equation.  In  another  tr< 

the  northwest  no  echo  was  observed  where  the  satellite  ESI  was  10.  he-.  * 

P 

vations  show  the  difficulties  of  comparing  ESI^  satellite  values  with  are..-  . t. 
echo  seen  only  in  a horizontal  plane. 


5.4.3  THE  14  AUGUST  197  6 CASE 

In  this  case,  the  cloud  sampled  was  an  active  shower  having  a top  at  ' 
rain  to  the  ocean  surface.  A racetrack  aircraft  ascent  was  made  along  'or  v 
(true)  radial  from  Kwajalein.  The  southern  end  of  the  track  was  inside  'hr 
The  northern  end  of  the  track,  during  ascent,  was  often  outside  the  cloud  >il  -i  • ■ 
in  the  altitude  range  of  9 and  13  km  where  cloud  extended  north  of  the  flig^*  • . - 
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Particle  classification  was  made  from  the  PMS-21)  samples.  I’MS-ll)  mass 
estimates  were  then  computed  and  plotted  as  shown  in  Figure  HO.  Observer  notes 
are  appended.  The  solid  profile  line  is  intended  to  represent  the  mean  mass  over 
a N-S  distance  of  about  22  km  or  three  satellite  grid  squares.  An  I0SI^  of  11.  6 was 
calculated  from  this  profile. 


Figure  HO.  Mass  Profile  and 
Observer  Notes  Through 
Active  Cb  Near  Kwajalein, 

2 103  to  2 1427.,  14  August 
1076.  Points  are  PMS-1D 
system  measurements.  Points 
with  crosses  (x!  are  averages 
over  more  than  one  sampling 
run 


Unfortunately,  the  NOAA  4 satellite  data  could  not  be  positioned  relative  to  the 
aircraft  with  certainty.  The  Kwajalein  WSR-57  radar  data  were  used  to  make  the 
adjustment.  Isolated  echoes  to  the  SSF  and  NNK  of  the  radar  were  matched  to  ap- 
parent cloud  centers  in  the  satellite  data.  The  resulting  best  fit  of  the  satellite 
data  and  radar  patterns  is  shown  in  Figure  HI.  Radar  pictures  at  16 -min  intervals 
showed  that  the  echoes  were  moving  from  162"  at  H3  km  hr  In  spite  of  careful 
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pattern  interpolation  during  the  sounding  period,  2104  to  2142/.,  few  echoes  were 
found  where  the  aircraft  was  supposed  to  be  sampling.  Some  uncertainty  remains 
since  larger  echoes  are  expected  from  such  a vigorous  cloud  system. 


Figure  31.  ESIp  (in  tenths),  WSK-57  Radar  Echoes  at  213  1/  (stippled)  and 
Aircraft  Sounding  Area  (heavy  box).  Satellite  data  at  2 138/  14  August  107  6 


Average  ESI^  over  the  three  satellite  squares  is  5.4  or  about  half  the  value 
found  by  aircraft.  On  the  other  hand,  if  the  aircraft  sounding  was  located  in  the 
same  size  area  but  one  square  to  the  left  and  two  squares  down  an  average  of  7. 8 
is  found. 

1 

5.5  Effect  of  Slant  Viewing  from  a Satellite 

Reference  was  made  in  Section  4.  6 to  the  utility  of  the  OOES-West  satellite 
data  for  delineating  areas  of  significant  ESI^  in  the  Kwajaletn  area.  The  OOES-West 
satellite  views  Kwajalein  from  an  angle  of  65°  off  vertical.  Hence,  the  atmospheric 
path  from  Kwajalein  to  OOES-West  is  more  than  twice  as  great  as  the  path  to  a 
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satellite  directly  overhead  and  atmospheric  attenuation  of  reflected  sunlight  or 
terrestrial  emission  could  be  significantly  greater  for  the  longer  paths.  Tests  were 
conducted  to  insure  that  the  ESI^  prediction  equations  were  useful  for  high  viewing 
angles.  An  example  is  shown  here  to  illustrate  the  effects  of  extreme  slant  viewing. 

In  order  to  simulate  conditions  when  the  Kwajalein  area  is  viewed  for  planning 
operations,  a cloud  system  was  viewed  from  GOES -East  looking  westward  at  about 
the  same  angle  that  GOES-West  looks  when  viewing  Kwajalein.  A time  was  chosen 
when  the  solar  zenith  angle  from  the  cloud  was  also  about  the  same  as  during  opera- 
tional usage.  The  angles  are  defined  in  Figure  4.  The  view  from  GOES -East  was 
then  compared  with  the  view  from  GOES-West.  In  this  test  the  zenith  angle  from 
the  cloud  to  the  eastern  satellite  was  56"  as  opposed  to  65°  from  Kwajalein  to  the 
western  satellite.  The  zenith  angle  from  the  same  cloud  system  to  the  western 
satellite  was  17°.  The  view  from  the  eastern  satellite  was  at  1600Z  and  the  view 
from  the  western  satellite  was  at  1615/..  In  this  test,  the  pattern  of  ESI^  computed 

from  the  western  satellite  is  considered  ESI  truth  since  it  had  viewing  conditions 
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close  to  those  in  the  development  sample.  It  is  shown  in  Figure  32.  The  ESI^ 
pattern  as  seen  slantwise  from  the  eastern  satellite  is  shown  in  Figure  33. 

Most  obvious  in  comparing  Figures  32  and  33  is  the  difference  in  E-W  scale. 
Secondly,  clouds  viewed  slantwise  are  displaced  westward  a distance  proportional 
to  their  height.  Clouds  in  the  vicinity  of  8.5N,  124.  OW  in  Figure  32  were  about 
13  km  high.  If  the  top  western  edge  was  sharp,  simple  geometry  places  the  edge 
20  km  further  west  than  its  true  position  in  the  slant  view.  Observed  displacement 
is  about  0.4"  or  44  km.  The  discrepancy  may  be  due  to  not  knowing  the  exact  nature 
of  the  cloud  edge,  that  is,  how  far  did  the  anvil  extend  westward  and  how  diffuse  was 
the  edge. 

In  Figure  32.  a large  area  with  ESI^  less  than  1.  0 is  northwest  of  a line  between 
8.  ON,  122.  3Wj  and  7 . 3 N,  and  123.  OW.  In  the  slant  view.  Figure  33,  the  same  area 
tends  to  have  ESI^  slightly  above  1.  0 with  many  fewer  zeros.  On  the  other  hand, 
areas  of  high  values  are  similar  from  both  angles.  For  example,  the  area  greater 
than  4.0,  centered  near  8.  5N,  123.  8W  in  Figure  32  is  about  the  same  in  size  in  the 
slant  view  of  Figure  33.  Peak  values  are  6.  2 vertical  vs  5.  4 slantwise.  Near  0.  8N, 
122.  2W  peak  values  are  10.0  vertical  vs  7.7  slantwise. 

Of  course  a comparison  of  this  type  assumes  correct  calibrations  of  both  the 
visible  and  IR  on  both  satellites  and  proper  normalization  factors  for  the  sun  angles. 
The  calibrations  were  tested  by  simultaneous  view  of  the  same  cloud  located  near 
the  equator  and  mid-way  between  the  satellites  when  the  sun  was  at  local  noon. 

These  comparisons  suggest  that  the  calibrations  do  not  yield  strictly  compatible 
data.  GOES-EAST  IR  temperatures  appeared  to  be  about  5"  lower  than  GOES-West 
at  200“K  and  l"  lower  at  200"K.  Calibrated  brightness  of  GOES-East  averaged 
2 percent  higher  than  Goes -West  with  no  discernible  trend  between  high  and  low 
vali.  s. 
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In  summary,  the  slant  view  gave  ESIp  about  20  percent  too  low  where  a more 

or  less  uniform  field  of  high  (3  to  10)  ESIp  was  observed.  Where  the  cloud  elements 

were  small  and  did  not  fill  the  sensor  field  of  view  in  the  slant  view,  the  ESI  values 

P 

were  low.  This  condition  can  be  detected  in  the  slant  view  by  a few  isolated  small 
ESI  values.  Widespread  small  values  of  ESIp  in  the  slant  view  were  probably  a 
field  of  broken  clouds  or  layer  of  thin  clouds  and  the  ESIp  values  were  higher  than 
those  in  the  vertical  view. 
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Figure  33.  ESIp  (in  tenths)  Pattern  Calculated  from  GOES  E,  56°  Zenith  Angle  to 
Satellite.  1600Z  25  .July  1977 


5.6  Satellite  v*  AKGL-ll  Erosion  Estimates 

Satellite  Estimates  of  ESI  were  compared  to  ESI  estimated  from  the  AFGL  II 

P P 

model  discussed  in  Section  1.  2.  The  comparison  is  made  possible  by  the  fact  that 
AFGL  II  hydrometeor  profiles  were  provided  every  three  hours.  Therefore,  once 
per  day  satellite  estimates  over  the  11  AFGL  II  stations  could  be  compared  to  the 
AFGL  II  estimates  closest  in  time  and  always  be  within  1.  5 hr  of  AFGL  II  esti- 
mates. As  in  earlier  applications  over  Europe  and  Asia,  the  satellite  estimates 
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were  based  on  averages  of  5 X 5 arrays  of  data  points  representing  areas  of  about 
50  X 50  km.  The  arrays  were  centered  as  close  as  possible  to  the  AFGL  II  station 
locations.  The  two  climatologies  had  an  overlapping  period  of  April  1973  to  Jan- 
uary 1974  . May,  July,  October,  and  January  were  chosen  as  representative 
seasonal  months  and  comparisons  were  made  for  all  11  stations  for  all  the  days  of 
those  months.  The  maximum  number  of  pairs  to  compare  would  be  1364;  however, 
only  1113  were  found  due  to  missing  satellite  data  for  some  dates  and  locations.  In 
the  months  May,  July,  and  October,  satellite  estimates  of  ESIp  were  from  the 
TR-  + T?^  equation  while  in  January  they  were  from  the  equation. 

ESI  values  were  also  calculated  from  the  AFGL  II  hydrometeor  profiles  using 
P 

two  different  methods.  Two  methods  are  available  since  the  AFGL  II  model  pro- 
vides not  only  estimates  of  hydrometeor  mass  densities  but  also  estimates  of  the 
percent  horizontal  coverage  of  cloud  and  precipitation  layers.  In  the  first  method, 
hydrometeor  density  was  simply  multiplied  by  the  fraction  of  horizontal  coverage 

before  ESI  values  were  calculated.  These  ESI  values  represent  areal  averages 
P P 

and  can  be  compared  directly  to  the  satellite  estimates  which  are  also  areal  aver- 
ages. In  the  second  method,  ESIp  values  were  calculated  assuming  that  all  cloud 
layers  had  100  percent  coverage.  These  values  describe  the  worst  case  situation 
of  a reentry  vehicle  encountering  clouds  at  all  the  cloud  layers.  These  values  can 
be  compared  to  the  satellite  estimates  if  the  satellite  estimates  are  increased  by 
the  ratio  of  the  Method  2 ESIp  divided  by  the  Method  1 ESIp.  For  either  method, 

AFGL  II  ESI  values  before  and  after  the  satellite  pass  time  were  linearly  inter- 
P 

polated  to  the  time  of  the  satellite  pass  before  comparison  to  the  satellite  values. 

Exceedance  statistics  for  AFGL  II  and  satellite  ESIp  values  in  four  months 

combined  are  compared  in  Figures  34  and  35.  The  plotted  lines  give  the  percentages 

of  time  when  ESI  thresholds  are  exceeded.  The  percentage  of  exceedance  de- 
P 

creases  as  ESI  increases.  Figure  34  has  the  results  for  Method  1,  in  which  re- 
P 

duced  AFGL  II  values  are  compared  to  satellite  values.  AFGL  II  percentages  are 
greater  than  satellite  percentages  for  ESIp  values  of  2 and  greater.  Inspection  of 
plots  for  individual  months  revealed  a general  tendency  for  greater  AFGL  II  per- 
centages above  an  ESI  value  of  2 with  the  extreme  differences  in  the  month  of  July. 

A definitive  explanation  of  the  July  differences  is  lacking  although  two  possibilities 
clearly  exist.  First,  the  50  X 50  km  averages  of  satellite  data  may  not  properly 
represent  the  convective  weather  of  July  and  thereby  underestimate  ESIp  as  dis- 
cussed in  Section  5.  4.  1.  On  the  other  hand,  the  second  possibility  is  that  the 
AFGL  II  data  include  surface  reports  of  cloud  coverage  which  overestimate  the 
layer  coverage  of  convective  clouds  since  the  observer  estimates  of  sky  coverage 
tend  to  combine  bases  and  sides. 
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Figure  35  has  the  results  for  Method  2,  in  which  maximum  AFGL  II  values  of 
ESI  are  compared  to  increased  values  from  the  satellite.  The  percentages  are 
moderately  higher  than  the  corresponding  Method  1 percentages.  In  Figure  35, 
AFGL  II  percentages  are  remarkably  close  to  the  satellite  percentages  for  ESI^ 
thresholds  of  five  or  less.  For  six  and  above  AFGL  II  percentages  are  higher,  as 
they  were  in  Method  1. 

These  comparisons  are  considered  to  be  good,  taking  into  account  that  the 
satellite  and  AFGL  II  estimates  are  from  two  completely  different  approaches  each 


with  substantial  uncertainty. 
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Figure  35.  Method  2 Exceedance  Statistics  for  AFGL  II  and  Satellite  ESI  Values. 
Satellite  values  are  increased  according  to  Method  2 (see  text),  AFGL  II  p values 
are  unchanged 


5.7  Conclusions  on  Verification 

An  ideal  experiment  for  verification  of  the  erosion  estimates  from  satellites 
would  have  the  following  elements.  First,  one  or  more  aircraft  would  be  available 
for  rapid  ascents  and  descents  to  measure  particle  size  distributions  and  phase  as 
a function  of  altitude  and  sufficient  for  numerical  models  of  reentry  vehicle  erosion. 
Then,  in  order  to  bridge  the  gap  between  the  very  small  volumes  sensed  from  air- 
craft and  the  large  volumes  sensed  from  satellites,  active  systems  such  as  radars 
or  lidars  would  scan  the  volume  sensed  by  the  satellite.  Moreover,  the  satellite 
measurements  would  be  calibrated  to  precise  physical  standards  which  could  be 
updated  at  the  time  of  cloud  observations.  The  satellite  measurements  would  be 
precisely  navigated  and  colocated  with  the  other  observations.  Given  these  capabili- 
ties, simultaneous  observations  of  clouds  would  be  taken  in  a great  number  of  cases 
for  different  cloud  types,  thicknesses,  coverages,  rainfall  rates,  temperature  pro- 
files, sun  and  observation  angles,  etc.  Given  these  data,  relationships  between 
satellite  measurements  and  erosion  parameters  could  be  checked  with  confidence. 
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None  of  the  cases  in  this  report  meets  all  or  even  most  of  the  elements  of  the 
ideal  experiment.  The  chief  sources  of  uncertainty  varied  from  case  to  case  de- 
pending on  how  well  the  various  instruments  performed,  how  quickly  clouds  were 
advecting  through  the  volume  sampled  by  the  satellite,  and  other  factors.  In  most 
cases,  the  hydrometeor  profiles  were  likely  to  have  more  error  than  the  satellite 
measurements.  All  of  the  41  cases  discussed  in  the  report  have  hydrometeor  pro- 
files based  on  aircraft  measurements  without  the  benefit  of  quantitative  radar 
measurements  in  3 dimensions.  The  reason  for  the  lack  of  radar  measurements  is 
simple.  Whenever  suitable  weather  existed  over  a calibrated  radar,  the  range  was 
reserved  for  erosion  tests.  Even  if  the  ranges  were  available,  the  long  delays  in 
waiting  for  suitable  weather  over  the  range  at  satellite  pass  time  would  probably  not 
allow  a sufficient  number  of  cases  for  a statistical  approach.  On  the  other  hand, 
the  aircraft  has  the  advantage  over  radars  by  measuring  particle  sizes  and  phases, 
and  in  sampling  thin  layers  of  cirrus  which  return  undetectable  echoes  to  the  radar. 
The  technology  of  particle  measurement  by  aircraft  was  undergoing  many  refine- 
ments during  the  three  years  in  which  data  were  taken  under  satellites  so  that  the 
latest  hydrometeor  profiles  were  greatly  improved  over  the  earliest.  When  both 
ascents  and  descents  were  made  by  the  Lear  aircraft  in  large  mid-latitude  cloud 
systems,  ESI^  values  varied  by  only  20  to  30  percent,  which  is  an  encouraging 
number. 

The  satellite  instruments  appear  to  have  reasonably  small  errors  according 
18 

to  Conlan.  This  claim  is  not  surprising  since  the  satellite  instruments  were  an 
operational  series  which  was  not  undergoing  development  changes.  When  the  rms 
error  of  2.02°K  in  the  satellite  IR  temperature  is  introduced  in  Eq.  (1)  of  Table  6, 
the  corresponding  ESI^  error  is  only  0.33  or  6 percent  in  the  worst  of  41  cases. 
Similar  errors  were  found  converting  temperature  to  altitude  (Table  14).  Visible 
data  are  subject  to  greater  uncertainty  due  to  corrections  for  solar  illumination  and 
bidirectional  reflectance;  however,  this  uncertainty  is  minimized  by  the  fact  that 
visible  data  are  weighted  much  less  than  IR  data  in  the  erosion  prediction  equations. 
No  allowance  was  made  for  background  brightness  or  temperature,  which  is  a source 
of  uncertainty  for  thin  clouds  or  partial  cloud  coverage. 

The  errors  in  earth  location  of  satellite  data  introduce  uncertainty  which  may 
exceed  the  instrumental  error  for  some  cases.  Despite  procedures  to  adjust  the 
archive  mapping  some  uncertainty  remains  as  to  the  exact  boundary  of  satellite  data 
sr  Me  for  comparison  to  the  aircraft  profiles  and  leads  to  uncertainty  in  the  areal 
averages  of  the  satellite  data.  The  uncertainty  increases  when  satellite  data  are 
averaged  over  smaller  areas.  Compared  to  the  satellite,  aircraft  navigation  was 
excellent  and  generally  needed  no  adjustment. 

The  uncertainty  in  applying  prediction  equations  developed  over  70  X 70  km 
averages  of  satellite  data  in  temperate  latitudes  to  tropical  cloudiness  near 
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Kwajalein  was  tested  in  three  cases.  For  these,  the  maximum  resolution  data  of 
7.5  X 7. 5 km  from  the  NOAA  archive  were  used  to  better  characterize  the  fine- 
scale  cloudiness.  In  the  case  of  a very  small  central  storm  core  which  reached  the 
surface  the  satellite  underestimated  the  core  ESI^  by  six  times.  Estimates  through 
thin  portions  of  the  anvil  were  in  close  agreement.  In  the  case  of  a large  dying 
shower  cloud  whose  virga  only  reached  downward  to  about  6 km  MSL,  satellite 
estimates  were  close  to  the  aircraft  measurements.  A third  case  must  be  viewed 
with  uncertainty  due  to  navigational  problems.  This  cloud  was  deep  and  active  with 
rain  reaching  the  surface.  The  satellite  estimates  were  from  50  to  60  percent  of 
the  A/C  estimates.  These  tests  indicated  that  reasonable  values  of  ESIp  are  given 
by  the  prediction  equation  when  maximum  resolution  data  are  used  without  averaging. 

For  all  of  the  erosion  parameters.  ESI^,  AS,  ADI,  and  AR,  multiple  regression 
correlation  coefficients  in  the  range  of  0.7  to  0.  8 were  consistently  found  from  the 
sample  of  41  cases.  The  coefficients  imply  that  a highly  significant  50  to  65  percent 
of  the  sample  variance  is  explained  by  the  equations.  When  the  41  case  sample  was 
randomly  divided  into  two  samples,  quite  similar  prediction  equations  were  found 
for  each  sample.  Therefore,  we  do  not  expect  that  larger  samples  or  a repeat  of 
the  experiment  would  change  the  correlation  by  much.  We  do  not  claim  to  have 
climatologically  representative  numbers  of  cloud  types  in  our  sample,  that  is,  we 
are  probably  underrepresented  in  cirrus  cases  and  overrepresented  in  multilayer 
and  thick  cloud  systems.  On  the  other  hand,  samples  and  tests  have  been  made  over 
a great  range  of  cloud  conditions  observed  by  aircraft  and  infrared  and  visible 
measurements  by  satellites.  Finally,  the  basic  relationship  between  satellite 

measurements  and  erosion  estimates  has  held  up  remarkably  well  in  a variety  of  * 

tests.  Properly  calibrated  data  from  GOES  geostationary  satellites  has  been  sub-  J 

stituted  for  data  from  NOAA  polar  orbiting  satellites  with  no  noticeable  loss  of 

usefulness.  Satellite  estimates  of  erosion  have  compared  well  with  estimates  from 

such  varied  sources  as  aircraft  sampling,  radar  mapping,  and  human  interpretation 

of  surface  weather  reports  and  radiosonde  data. 

;J 


77 


! 


References 


1.  Smith,  R.  C.  , Capt  (1974)  Atmospheric  Moisture  Parameterization,  USAFETAC 

TN  74-1,  Washington,  D.  C.  19  pp. 

2.  Feddes,  R.G.,  Capt  (1974)  A Synoptic-scale  Model  for  Simulating  Condensed 

Atmospheric  Moisture.  USAFETAC  tN  74-1,  Washington,  D.  C.  , 21  pp. 

3.  Cunningham,  R.  M.  , and  Peirce,  R.  M.  (1974)  Environmental  Definition  Plan 

Cross  Section  Analysis,  AFGL  (unpub.  ms.  ). 

4.  Feteris,  P.  J. , Lisa,  A.  S. , and  Bussey,  A.  J.  (1975)  Environmental  Definition 

Program  Cross  Sectional  Analysis;  Summary  of  Data  and  Analysis  Tech 
niques.  AFCRL-TR-76-0602.  Final  Rept.  F19628-74-C-0073,  Environmental 
Res.  and  Tech.  , Inc.,  Concord,  Massachusetts  01742,  34  pp. 

5.  Gringorten,  1.  1.  (1976)  Areal  Coverage  Estimates  by  Stochastic  Modeling. 

AFGL-TR-76-0148,  ERP  No.  '573,  56  pp.  

6.  Liou,  K-N,  and  Stoffel,  T.  (1976)  Remote  Sensing  of  Cirrus  Cloud  Composi- 

tions From  Satellites,  AFGL-TR-76-0<)27,  81  pp.  University  of  Utah, 

Salt  Lake  City,  Utah  84112. 

7.  Feddes,  R.  G. , and  Liou,  K-N  (19  7 7)  Cloud  Composition  Determination  by 

Satellite  Sensing  Using  the  Nimbus  VI  High  Resolution  Infrared  Sounder, 
AFGL-TH-77-0123,  140  pp.  University  of  Utah,  Salt  Lake  City,  Utah  84 112. 

8.  Staelin,  D.  H. , Kunzi,  K.  F. , Pettyjohn,  R.  L.  , Poon,  R.  K.  L.  , and  Wilcox,  R.  W. 

(1976)  Remote  sensing  of  atmospheric  water  vapor  and  liquid  water  with  the 
Nimbus  5 microwave  spectrometer,  J.  Appl.  Meteor.  15:1204-1214. 

9.  Wilheit,  T.T.,  Chang,  A.  T.  C. , Rao,  M.S.V.,  Rodgers,  E.  B. , and  Theon,  J.S. 

(1977)  A satellite  technique  for  quantitatively  mapping  rainfall  rates  over 
the  oceans,  J.  Appl.  Meteor.  16:551-560. 

f 1 vw 

10.  Booth,  A.  L.  (1973)  Objective  cloud  type  classification  using  visual  and  infrared 

satellite  data.  Proceedings  of  the  3rd  Conf.  on  Probability  and  Statistics  in 
Atmospheric  Science,  pp  220-227. 

11.  Shenk,  W.E.  , Holub,  R.  J.  , and  Neff,  R.  A.  (1976)  A multispectral  cloud  type 

identification  method  developed  for  tropical  ocean  areas  with  Nimbus-3 
MRIR  measurements.  Mo.  Wea.  Rev.  104:284-291. 


References 


12. 

13. 

14. 

I 15. 


16. 


17. 


18. 

19. 

20. 

21. 

22. 


23. 

24. 

25. 

26. 
27. 


28. 


Reynolds.  IX  . and  Vonder  Haar,  T.  H.  (1973)  A comparison  of  radar-deter- 
mined cloud  height  and  reflected  solar  radiance  measured  from  the  geo- 
synchronous satellite  ATS-3,  J.  Appl.  Meteor.  E2: 1082- 1085. 

Griffith.  C.G. , and  Woodley,  W.  L.  (1973)  On  the  variation  with  height  of  the 
top  brightness  of  precipitating  convective  clouds,  J.  Appl.  Meteor. 
12:1086-1089. 

Kaveney,  W.  J. , Feddes,  K.  , and  Liou,  K-N  (197  6)  Statistical  Inference  of 
Cloud  Thickness  from  NOAA  4 Scanning  Radiometer  Data,  Al',GL-,fll-76- 
0 1 05 , 76  pp.  University  of  Utah,  Salt  Lake  City,  Utah  84 112. 

Martin,  D.  W. , and  Scherer,  W.  D.  (1973)  Review  of  satellite  rainfall  estima- 
tion methods.  Bull.  Amer.  Meteor.  Soc.  54:661-674. 


Griffith,  C.  G. , Woodley,  W.  L. , Grube.  P.  G. , Martin,  D.  W.  , and  Stout,  J. 
(1977)  Rain  Estimation  from  Geosynchronous  Satellite  Imagery,  Visible 
and  Infrared  Studies,  submitted  to  Monthly  Weather  Review  from  NHEML, 
NOAA,  Coral  Gables,  Florida. 

Bunting,  J.T.  (1976)  Cloud  Properties  from  Satellite  Infrared  and  Visible 
Measurements,  AFGI.-TR-76-0287,  6 pp. 

Conlan,  E.  F.  (1973)  Operational  Products  from  ITOS  Scanning  Radiometer  Data, 
NOAA  Technical  Memorandum  NESS  52,  57  pp. 

Sikula,  G.J.,  and  Vonder  Haar,  T.  H.  (1972)  Very  Short  Range  Local  Area 

Weather  Forecasting  Using  Measurements  from  Geosynchronous  Meteorologi- 
cal' S'iteiritTSr^Tnar  Report  For  Contract  F 1 5628-7 1 -£-0073,  ATCHL-76- 
DM 

Raschke,  E. , Vonder  Haar,  T.  11. , Pasternak,  M.  , and  Bandeen,  W.  R.  (1973) 
The  Radiation  Balance  of  the  Earth- Atmosphere  System  from  Nimbus  3 
Radiation  Measurements,  NASA  TN  D-724S),  71  pp. 

Cozzens,  D.  E.  (1974)  Technique  to  Earth  Locate  and  Analyze  Data  from 
Satellite  Sensors,  pp  13-36.  Burke,  L.  (1974)  Met  and  Math.  Anal.  Final 
Rept.  Regis  College.  AFCRL-TR-74-0558. 

Bristor,  C.  L.  (1975)  Central  Processing  and  Analysis  of  Geostationary  Satellite 
Data.  NOAA  Tech  Mem.  NESS  64,  Washington,  D.  C.  , 155  pp. 

Church,  J.  F. , Lt  Col,  Poes,  K.  K.  , and  Spatola,  A.  A.  (1975)  The  Continuous 
Aluminum-Foil  Hydrometeor  Sampler;  Design,  Operation,  Data  Analysis 
Procedures,  and  Operating  Instructions.  AFCRL-TR-75-0370,  70  pp. 

Plank,  V.G.  (1974)  Hydrometeor  Parameters  Determined  from  the  Radar  Data 
of  the  SAMS  Rain  Erosion  Program,  SAMS  Rept.  No.  2,  AFCRL-TR-74- 
0249,  86  pp. 

Hobbs,  P.  V. , Radke,  L.  F.  , and  Atkinson,  D.G.  (1975)  Airborne  Measurements 
and  Observations  in  Cirrus  Clouds,  AFCRL-TR-75-0249,  Sci.  Rept.  No.  1, 
University  of  Washington,  Seattle,  Washington,  117  pp. 

Berthel,  R.  O.  (1975)  A Climatology  of  Selected  Storms  for  Wallops  Island, 
Virginia.  1971-1975,  SAMS  Rept.  No.  4.  AFGL-TR-76-0118.  82  pp. 

Booker,  R.  D.  , and  VVindes,  J.  (1977)  HA  WADS  Equipment  Description  and 
Operational  Manual,  Aeromet,  Inc.  , At'GL-TR-77-0066,  115  pp. 

Knollenberg,  R.G.  (1975)  The  Response  of  Optical  Array  Spectrometers  to 
Ice  and  Snow:  A Study  of  f* robe  Size  to  Crystal  Mass  Relationships, 

Sci.  ftept.  No.  1,  Particle  Measuring  Systems,  Inc.,  Boulder,  Colo rado 
80301,  AFCRL-TR-75-0494,  70  pp. 


1 


80 


References 


Knollenberg,  R.  G.  (1976)  The  Response  ofOptical  Array  Spectrometers  to  Ice 
and  Snow;  a Study  of  2-D  Probe  Area-to-Mass  Relationships,  Final  Rept. 
Particle  Measuring  Systems,  Inc.,  Houlder,  Colorado  8(1301,  AFGL-TR- 
76-0273. 


30.  Serbagi,  R.C.,  Rodenhiser,  P.  M.  , and  Kaplan,  F.  B.  (197  4)  Warm  Fog  Dissi- 
pation  Models  and  Control  and  Processing  Programs  for  Airborne  Data 


bystem,  AFCKL-T  K -74  -04  I X,  164  pp,  Digital  Frog.  Services,  Inc., 
Waltham,  Massachusetts  02154. 

31.  Belsky,  L.  E.  , Kaplan,  F.  B.  , and  Rodenheiser,  P.  M.  (197  5)  Development  and 
Application  of  Data  Processing  Techniques  and  Analytic  Procedure  to 


. pp,  Digital  Prog.  Services, 

Inc. , Waltham,  Massachusetts  02154. 

32.  Belsky,  L.  F. , Francis,  M.W.,  Kaplan,  F.  B.  , and  O'Neil,  J.  E.  (1976)  Con- 

tinuation  of  Development  and  Application  of  Data  Processing  Techniques  and 
Analytic  Procedures  to  Cloud  Physics  Data,  AFGL -Tft -1  6-0 182,  198  pp. 
Digital  Prog.  Services,  Inc. , Waltham,  Massachusetts  02154. 

33.  Spiegler,  D.B.,  and  Fowler,  M.  G.  (1972)  Four  Dimensional  World  Wide 

Atmospheric  Models  (surface -25  km  altitude)  NASA  CR-2082,  Environmental 
Res.  and  Tech.  , Inc.,  Concord,  Massachusetts  01742,  65  pp. 

34.  Crutcher,  H.  L. , and  Meserve,  J.  M.  (1970)  Selected  Level  Heights,  Tempera- 

ture  and  Dew  Points  for  the  Northern  Hemisphere,  NaVaIR  50-1C-52. 


Append'*  A 

Computer  Processing  of  Archived  Satet'ite  Oeu  • 
Climatotogicei  Appt.tr’  ' 

The  production  of  erosion  climatologies  over  target  areas  and  test  • , • . 

required  extensive  computer  programming  and  execution  time.  The-.-  < . 

of  programs  was  used  for  each  climatology.  The  only  difference  was  the  • . if. 
tion  of  a different  box  of  rows  and  columns  within  the  NOAA  Norther-  H«- 
satellite  archive.  For  the  Eurasian  climatology,  the  box  was  large  »r- 
data  points  covering  most  of  Europe  and  Asia.  With  arrays  of  this  sire.  • , 
convenient  to  work  with  only  4 days  at  one  time.  For  the  test  range  • . ... 

the  boxes  were  much  smaller  with  20  X 25  data  points  covering  impact  . ,• 

Kwajalein  Atoll  and  Wallops  Island,  Virginia.  With  these  smaller  :irru-  ■ . . 
convenient  to  work  with  8 days  and  both  areas  at  one  time. 

The  following  list  gives  the  titles  and  functions  for  programs  to  nr-.  *•-- 
archived  NOAA  satellite  data  to  produce  erosion  climatologies.  The  \t.  \ 
was  described  in  Section  2.  1.  An  estimate  of  CP  time  on  the  A F < .1  IH 
puter  for  the  Eurasian  climatology  is  given  for  the  programs. 

Program  Function 

CORFACT  Determines  orbits  and  calculates  viewing  i-.g  •• 

each  10  X 10  array  of  data  for  each  day  on  t h*  \<  \ . n 

polar  sterographic  archive  tapes.  Selects  .:••  • 
correction  factors  for  visible  data  from  tab >-  •• 
bidirectional  reflectance.  (380  sec  dav 
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Procram 


K unction 


SR PACK 


SRCORR 


SHKSI 


SR  ESI  P 


CO/.I’T 


PRORMP 


rsr^r'- 


Pnpacks  data  from  SR  archive  tapes  and  saves  data 
in  predetermined  boxes.  Checks  for  bad  records. 

(350  sec /day) 

Merges  data  from  2 tapes  to  eliminate  the  "dateline" 
over  Europe  in  the  SR  archive  tape3  and  multiplies 
all  visible  data  by  appropriate  correction  factors. 

(500  see/dav) 

Averages  visible  and  IR  data  over  5X5  arrays. 

Tests  for  bad  or  missing  data.  Output  options  are 
erosion  parameter  estimates.  (170  sec/day) 

Calculates  estimates  of  WSI  . ADI.  AS,  AH,  and 

P 

also  total  EWC  from  average  visible  and  IR  for 
summer  months  and  IR  and  temperature  profiles 
for  winter  months.  (30  sec/day) 

Generates  computer  plots  of  areas  where  erosion 
parameters  exceed  certain  thresholds.  National 
boundaries,  coastlines,  and  a latitude-longitude 
grid  are  plotted  for  reference.  See  Figure  2 of  Annex, 
(30  sec/day) 

Computes  probabilities  of  exceeding  erosion 
parameters  at  all  locations  for  a given  month. 

The  probabilities  are  used  to  adjust  lor  missing 
lata.  (5  tec  I day) 
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